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We describe the design, construction and performance of the PI- 
BETA detector built for the precise measurement of the branching 
ratio of pion beta decay, tt"*" —>■ ir^e'^Ue, at the Paul Scherrer In- 
stitute. The central part of the detector is a 240-module spherical 
pure Csl calorimeter covering ~37rsr solid angle. The calorimeter 
is supplemented with an active collimator /beam degrader system, 
an active segmented plastic target, a pair of low-mass cylindrical 
wire chambers and a 20-element cylindrical plastic scintillator ho- 
doscope. The whole detector system is housed inside a temperature- 
controlled lead brick enclosure which in turn is lined with cosmic 
muon plastic veto counters. Commissioning and calibration data 
were taken during two three-month beam periods in 1999/2000 
with vr"*" stopping rates between 1.3-10^ tt~^ /s and l.STO^ 7r~^ /s. We 
examine the timing, energy and angular detector resolution for pho- 
tons, positrons and protons in the energy range of 5-150 MeV, as 
well as the response of the detector to cosmic muons. We illustrate 
the detector signatures for the assorted rare pion and muon decays 
and their associated backgrounds. 

PACS Numbers: 29.40.Gx, 29.40Vj, 29.85.+C 
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tillators; Tracking and particle identification; Detector design, con- 
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1 Introduction 



The PIBETA collaboration [1] has conducted a program of measurements at 
the Paul Scherrer Institute (PSI) with the ultimate goal of determining pre- 
cisely the branching ratio of pion beta decay (tt/?), R{'k^ n^e^Ue). tt/? decay 
is one of the most fundamental semileptonic weak interactions. It is a tran- 
sition between two spin-zero members of an isospin triplet, and is therefore 
analogous to super-allowed pure Fermi transitions in nuclear beta decay [2]. 
Due to its relative simplicity and, consequently, its accurate theoretical de- 
scription, Fermi beta decay stands today among the prominent successes of 
the Standard Model of electroweak interactions [3] . 

The most recent and most precise determination of the 7r/3 decay rate I/ttt/j is 
due to McFarlane and collaborators [4] , who performed a decay-in-flight mea- 
surement using the LAMPF spectrometer. Their result, 0.394±0.015 s~^, is 
in good agreement with the Standard Model prediction of 0.4027±0.0018 s~^ [5], 
but with a ~ 3.8 % uncertainty it does not test the full extent of the radiative 
corrections which stand at ~ 3 % [6,7]. When combined with the tt"*" lifetime [8], 
the rate measured by McFarlane et al. gives the 7r/3 branching ratio relative 
to aU pion decays of Kp/Ttotax = (1-026 ± 0.039) • 10"^ 

In the first phase of the PIBETA experiment we have collected ~60,000 nf] 
events, and we expect to extract the branching ratio with ~0.4% statistical 
and comparable systematic uncertainty. We have cross-calibrated the abso- 
lute TT"*" beam intensity using the measured rate of the tt"*" — > e'^u decays 
which is known with the combined statistical and systematic uncertainty of 
~0.40% [9,10]. The proposal for a second phase of the experiment, intended 
to reduce the overall uncertainty of both the tt/S and the n ^ eu measurements 
to ~0.2% is now under consideration. 

The PIBETA apparatus is a large solid angle nonmagnetic detector optimized 
for measurements of photons and electrons in the energy range of 5-150 MeV. 
The main sensitive components of the apparatus, shown and labelled in Fig. 1, 
are: 

(1) a passive lead collimator, PC, a thin forward beam counter, BC, two 
cylindrical active collimators, ACi and AC2, and an active degrader, AD, 
all made of plastic scintillators and used for the beam definition; 

(2) a segmented active plastic scintillator target, AT, used to stop the beam 
particles and sample lateral beam profiles; 

(3) two concentric low-mass cylindrical multi-wire proportional chambers, 
MWPCi and MWPC2 surrounding the active target, used for charged par- 
ticle tracking; 

(4) a segmented fast plastic scintillator hodoscope, PV, surrounding the 
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MWPCs used for particle identification; 

(5) a high-resolution segmented fast shower Csl calorimeter surroimding the 
target region and tracking detectors in a near-spherical geometry; 

(6) a cosmic muon plastic scintillator veto counters, CV, around the entire 
apparatus (not shown). 

All detector components listed above, together with the photomultiplier tube 
(PMT) delay cables, high voltage (HV) supplies, MWPC instrumentation and 
gas system, fast trigger and digitizing electronics, two front-end computers, 
and temperature control system are mounted on a single platform and can 
be moved as a self-contained unit into the experimental area (Fig. 2). When 
the detector platform is precisely positioned and surveyed with respect to 
the beam line, electrical power and Ethernet connections make the detector 
immediately operational. 

The structure of this article is as follows. The secondary beam line elements, 
the beam tunes and the overall layout of the experimental area are presented 
in Sec. 2. In Sees. 3-8, we describe the design, fabrication and quality control 
of the individual detector components enumerated above. In particular, the 
Csl calorimeter is covered in detail in Sec. 6. Sections 9-12 cover the PMT 
waveform digitizing system, the fast trigger electronics, and the data acquisi- 
tion system with user analysis software. The calibration procedures and the 
performance of the detector are demonstrated in Sees. 13-16. There we discuss 
detector timing, energy and angular resolution, particle identification, track 
reconstruction, and signal-to-background ratios for the selected pion and muon 
decay processes. Finally, a short overview of the Monte Carlo (MC) detector 
simulation coded using the GEANT package [11] is given in Sec. 17. At the end, 
in Sec. 18, we review the radiation resistance of the active detector elements. 



2 Beam line and experimental cirea 



The experiment was mounted at the Paul Scherrer Institute (PSI) in ViUigen, 
Switzerland. The isochronous separated sector cyclotron at PSI accelerates 
protons to an energy of 590 MeV. The 1.5 mA proton beam is transported 
along the primary proton channel to two target stations where pions and 
muons are generated and transported via secondary beam-lines to the exper- 
imental areas. The accelerator operates at the frequency of 50.63 MHz pro- 
ducing a microscopic beam structure of 1 ns wide proton pulses separated by 
19.75 ns. 

The PIBETA apparatus is set up in the PSI ttEI experimental area whose 
16 m long beam line is designed to supply intense low energy pion beams 
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with good momentum resolution [12,13]. Pions are extracted from the second 
target station E at an angle of 8° with respect to the incident protons. During 
the PIBETA measurements, this graphite target was 60 mm in length along 
the axis of the proton beam. Operating in a high-flux optical mode, the ttEI 
beam line can deliver a pion beam with a maximum momentum of 280 MeV/c, 
a Pull- Width-Half-Maximum (PWHM) momentum resolution of < 2% and an 
accepted production solid angle of 32 msr. The primary proton current in the 
ring cyclotron during the PIBETA data acquisition periods in years 1999-2001 
was 1.6 mA DC on average. 

Fig. 3 shows a sketch of the fixed beam line elements between the pion produc- 
tion target and the last pair of focusing magnets (which are embedded in the 
concrete shielding of the proton channel) . The channel incorporates three slit 
systems that are used to control the beam intensity either by coUimating the 
vertical size of the beam (FSH51) and the beam halo (FS51), or by reducing 
the momentum band acceptance (horizontal jaws FSH52), hence improving 
momentum resolution of the transported beam [14]. The bending dipole mag- 
net ASZ51 located in the front section of the beam line is used to select the 
beam momentum. 

We have developed a 113.4 MeV/c tt"*" beam tune with FWHM momentum 
resolution Ap/p < 1.3% and maximum nominal 7r+ beam intensity /,r = 
1.4- 10*^ 7r+/s, reached at the full cyclotron current of 1.7 mA. Beam tunes with 
TT"*" fluxes below 200k /s have a narrower momentum acceptance of Ap/p < 
0.28% (FWHM). The choice of a particular beam momentum is governed 
by the need for good time-of-flight (TOP) separation of pions, positrons and 
muons between the production target E and the first beam defining counter 
BC, as well as between the beam counter BC and the stopping target AT (see 
below). 

To reduce positron contamination even further, a 4 mm thick carbon de- 
grader is inserted in the middle of the ASY51 dipole magnet. The momentum- 
analyzed pions and positrons have different energy losses in the carbon ab- 
sorber, and are therefore spatially separated in a horizontal plane. Unfor- 
tunately, this also broadens the beam phase space. We have used TRANS- 
PORT [15] and TURTLE [16] beam transport codes to develop a nontradi- 
tional beam optics with foci in both the horizontal and vertical planes at the 
PSH52 momentum-limiting slit. The resulting beam tune reduces the phase 
space broadening introduced by the carbon degrader. A significantly higher 
luminosity at the PIBETA target position is thus achieved and the pions are 
stopped in a laterally smaller region. 

A simplified layout of the PIBETA apparatus is depicted in Fig. 4. A lead brick 
collimator PC with a 7 mm pin-hole located 3.985 m upstream of the detector 
center restricts the spatial spread of the incident tt"^ beam. The beam particles 
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are first registered in a 3 mm thick plastic scintillator BC placed immediately 
upstream of the collimator. The pions are subsequently slowed in a 40 mm 
long active plastic degrader AD and stopped in an active plastic target AT 
positioned at the center of the PIBETA detector. Fig. 5 shows the TURTLE- 
calculated momentum spectrum of the tt"*" beam incident on the front face of 
the degrader counter. The predicted vr"*" FWHM momentum spread at this 
point is 1.2MeV/c/113.4MeV/c~l.l%. 

The PSI surveying group measured the detector location using a transit theodo- 
lite. The upstream and downstream detector sides are levelled, with the center 
1506.0 ± 0.3 mm above the floor, corresponding to the nominal beam height 
in the vrEl channel. The measurement of the detector center is reproducible 
and the alignment of the target axis along the detector longitudinal axis is 
confirmed with an uncertainty of 0.3 mm. 

We used the OPTIMA control program [17] to adjust the currents in the 
dipole and quadrupole beam line magnets that steer and focus the tt"*" beam 
into the target. The goal was to achieve the smallest, most symmetric beam 
spot consistent with the high 7r+ intensity. The OPTIMA program allows a 
user to maximize an arbitrary experimental rate normalized to the primary 
cyclotron current by iteratively changing the settings of the magnetic elements. 
We chose to maximize the rate of four-fold coincidences between the forward 
beam counter BC, the degrader counter AD, the active target AT and the 
accelerator rf signal. These four signals are combined in a coincidence unit 
in such a way that their overlap signals a vr"*" particle stopping in the active 
target. Fig. 6 shows the relative timing of the vr^-stop trigger signals. The 
final optimized settings of beam line magnets are written in a disk file as well 
as saved and updated continuously in the online experiment database. 

The beam line shutter that regulates the beam's entry into the experimental 
area is computer operated. This option is very useful for restarting automatic 
data taking in instances when the beam goes down temporarily and detector 
operation has to be restarted from a remote location. 



3 Beam-defining detectors 



3. 1 Forward passive beam collimator 



The lead collimator PC is positioned immediately upstream of the fixed beam 
line's exit window, approximately at the focal plane between the last quadrupole 
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doublet (QSL53/54) and the final beam focusing triplet (QSK51, QSL55, 
QSK52). Only the beam pions pass through the collimator, while positrons 
(or electrons) stop in the collimator material, since at this point they are sep- 
arated from pions by ~ 40 mm in a horizontal plane by the bending magnet 
ASL51. 

This passive collimator consists of two stacked lead brick blocks with individ- 
ual dimensions of 250 x 250 x 50 mm^. Both pieces have a central hole with 
a diameter of 50 mm and a step bore extending the hole to 70 mm. Remov- 
able lead plugs with different apertures can be inserted into these openings. 
The upstream collimator plug used in production running has a central hole 
diameter of 7 mm. The hole in the downstream plug is 10 mm in diameter. 
The two collimator pieces are screwed together and then mounted as a unit 
onto a vacuum blind fiange which in turn couples to a beam window fiange. 
Therefore, the coUimator assembly always has a fixed position with respect 
to the beam vacuum tube. The complete collimator setup can be moved con- 
tinuously in both the horizontal and vertical directions with a quasi-rack and 
pinion mechanism in order to adjust its pin hole position relative to the pion 
beam axis. 

The pion beam transmittance is checked periodically by taking Polaroid pho- 
tographs of the beam spot at the vertical plane in front of the collimator's 
central hole. Fig. 7 shows a photograph from our experimental logbook. A one 
minute beam exposure reveals the resolved tt"*" and e+ beam spots, separated 
at this point by ~ 40mm. The circle centered on a pair of cross-hairs represents 
the position and size of the collimator opening. 



3.2 Forward beam counter 



The forward beam counter BC is the first active detector counter. It is placed 
immediately downstream of the lead collimator and just upstream of the verti- 
cal bending magnet SSL51. This counter tags beam particles that pass through 
the collimator. 

The central part of the beam counter is a rectangular piece of BICRON BC- 
400 plastic scintillator with dimensions 25 x 25 x 2mm^. The scintillator is 
optically coupled on all four edges to tapered acrylic light guides. One light 
guide is glued to the scintillator edge surface with BICRON BC-600 optical 
cement. The other three light guides are coupled to the scintillator via air 
gaps. Both the scintillator and the light guides are mounted inside a light- 
tight enclosure which consists of an aluminum frame with outer dimensions 
150 X 150x50 mm^, covered by two thin (30 /xm) aluminum windows. The frame 
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is attached to the lead colhmator and fixes the counter position. Mechanical 
feedthroughs at the four sides of the box hold the light guides and the magnetic 
shield cylinders of the PMTs, as seen in Fig. 8. 

Each of the four hght guides is air-gap coupled to a Hamamatsu R7400U mini- 
PMT. The PMT voltage dividers were custom designed at the University of 
Virginia to operate at typical counter rates > 2 MHz. The four analog PMT 
signals are electronically summed in a NIM LeCroy Model 428F linear fan-in 
unit. The two identical outputs representing summed analog BC signals are 
distributed to the fast trigger electronics and the digitizing electronics branch 
(see Sec. 9). Fig. 9 shows a pulse charge spectrum of BC signals digitized 
in a 10-bit CAMAC ADC LeCroy Model 2249 with a 25 ns integration gate. 
Accounting for the pion energy straggling in the thin counter, the 12.7% rms 
resolution of the ADC spectrum corresponds to a scintillator hght output of 
~ 103 photoelectrons/MeV. 



3.3 Active beam collimator pair 



Beam-related particles can reach the PIBETA calorimeter accidentally in co- 
incidence with a TT-stop signal. The accidental coincidences arise from: 

(1) beam halo particles entering the detector off axis; 

(2) beam particles scattering off detector support striictures; 

(3) beam pions undergoing hadronic interactions in air or in passive and 
active detector elements; 

(4) particles from pion and muon decays in flight occurring outside the active 
target volume. 

Two active beam collimators ACi and AC2 suppress backgrounds caused by 
detector hits that are not associated with a 7r-stop event. Both collimators are 
cylindrical rings made of 25.4 mm thick BICRON BC-400 plastic scintillator. 
The first (upstream) counter has an outer diameter of 120 mm and a central 
hole diameter of 50 mm. The second (downstream) active collimator has an 
outer diameter of 172 mm and a central hole diameter of 90 mm. The collimator 
ACi is placed around the beam vacuum tube, about 550 mm upstream from 
the center of the detector. The colhmator AC2 is attached to the aluminum 
ring mounting of the active degrader and covers its struts. It is located 310 mm 
upstream of the center of the detector. The collimator dimensions were chosen 
so as to ensure that beam tube elements and mechanical support structures 
are shielded without intersecting the calculated envelope of the beam. 
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The collimator rings arc flattened on the outer envelope at four points by 2 mm 
deep cuts, every 90 degrees. Light is collected by Hamamatsu R7400U mini- 
PMTs glued at each of these flat surfaces. The analog pulses from the four 
PMTs are summed in a NIM LeCroy Model 428F linear fan-in unit located in 
the electronics hut. This provides one common analog signal for each active 
collimator counter. 

Typical collimator counting rates at ~ 1 MHz stopping tt^ flux are 10 kHz and 
20 kHz for ACi and AC2, respectively. These comparatively low hit rates allow 
us to use the collimator signals for veto cuts in the offline analysis without 
introducing an additional dead time inefficiency. 



3.4 Active degrader 



In order to maximize pion beam transmission along the length of the beam 
line, pions reaching the central detector area have a relatively high momen- 
tum of 113.4MeV/c. They are subsequently moderated by the active de- 
grader AD. The degrader counter reduces the average pion kinetic energy 
from 40.3 MeV to 27.6 MeV. We continually monitor the ratio of counting 
rates t = A^(BC ■ AD)/A^(BC) to keep the beam transmittance stable. The 
theoretical value calculated for perfect beam transport with pions decaying 
in flight in a 3.82 m long region is 46.0 %. The measured ratio is in the range 
0.45 < t < 0.48, reflecting an efficient beam transport with minimal losses. 

The active degrader counter is made of BICRON BC-400 plastic scintillator 
and has the shape of a truncated cone. Fig. 10 shows the cross section of 
the degrader counter. This geometry ensures that the degrader 's downstream 
projection covers the whole target area (40 mm diameter) while particles en- 
tering parallel to the beam axis always traverse the same 30 mm scintillator 
thickness. Furthermore, the degrader's radial dimensions are compatible with 
the inner MWPC cylinder, which has a 90 mm minimal inner diameter. The 
outer surfaces of the degrader counter are therefore slanted and connect to 
four acrylic light guides, each equipped with a Hamamatsu R7400U PMT. 
This geometry is required so that the PMTs are located outside of the inner 
region, while the degrader projection matches the target area. The slant an- 
gle is small, optimizing light collection. The flsh-tail shaped light guides are 
bent into quarter circles. The light guides couple to the whole outer area of 
the scintillator. Each light guide ends up in a square cross section 6 x 6mm^, 
which matches the round window of the PMT. The photomultiplier tubes are 
coupled with an air gap. 

The degrader counter is supported by a plastic stand which also holds the 
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PMTs and voltage divider bases. It is mounted at the ends of four steel 
threaded rods protruding from the outer support structure of the detector 
into its center from the upstream side, leaving a 5 mm gap to the active target 
face. The inner active collimator is mounted to the same mechanical support 
structure and covers the rods in the direction of the beam. 

Fig. 11 shows the pulse-height spectrum of the active degrader for the vr^-in- 
beam trigger (see Sec. 10.2) taken with a short, 25 ns gate in a 10-bit CAM AC 
ADC unit. The spectrum has an excellent rms energy resolution of 7.8 %, cor- 
responding to photoelectron statistics of ~ 270photoelectrons/MeV. It easily 
identifies the events with two or more beam particles, from either a single 
beam pulse or two adjacent beam pulses. 



3.5 Active PIBETA stopping target 



The active target is a cylindrical plastic scintillator counter with a 50 mm 
length and a 40 mm diameter. The counter is segmented into 9 elements, as 
shown in Fig. 12. These 9 pieces are optically isolated by wrapping them 
individually in aluminized Mylar foil. The segments are pressed together and 
the whole assembly is wrapped with black plastic tape. Each target element 
is coupled to a miniature (8 mm photo-cathode) Hamamatsu R7400U PMT 
via a 60 mm long, tapered acrylic light guide. Fig. 13 is a photograph of the 
partially assembled target counter. 

The analog signal from each of the target segments is divided into two branches 
by a custom-made passive splitter. One side is discriminated in a Phillips 
Scientific PS 7106 16-channel time-over-threshold discriminator packaged in a 
single width CAMAC module. This output is digitized in a FASTBUS time-to- 
digital converter (TDC) LeCroy Model 1877, as well as counted with a scaler 
unit (PhiUips PS 7132H) that is read out every 10 s. The second branch is 
connected to an 12-bit FASTBUS analog-to-digital converter (ADC) LeCroy 
Model 1882F, gated with a 100 ns event trigger gate. 

Fig. 14 shows the FASTBUS ADC spectrum in a single target detector for 
7r"''-in-beam trigger events. The ~ 10 % rms width of the 7r+ stopping peak 
at 27.6 MeV is consistent with ~ 8photoelectrons/MeV. PMT gains of the 
target counters are equalized using a 4.1 MeV tt ^ /lu stopping muon line as 
a reference. 
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3.6 Beam stopping distribution and beam composition 



We have extracted the counting rates of vr-stop events and beam positrons in 
the individual target segments from the gain-matched ADC spectra. Fig. 15 
shows the 2-dimensional shape of the tt"^ stopping distribution using the count- 
ing rates of individual segments superimposed on an outline of the target. 
The TT"*" beam spot measured by this method agrees with the independent 
back-tracking tomography analysis that uses the tt"*" — > e^Ug positron tracks 
recorded by the MWPCs. The beam spot is almost symmetric and is cen- 
tered on the target, with horizontal and vertical profile rms widths of 7.6 mm 
and 8.4 mm, respectively. The 2-dimensional x-y beam profile is represented in 
our analysis by two separated, rotated, and then modulated Gaussian distribu- 
tions. The details of our algorithm are described in a separate publication [18]. 

The longitudinal distribution of stopping pions in the active target scintillator 
is calculated in a GEANT simulation and agrees very well with the back-tracking 
tomography reconstruction of the z coordinate beam profile. Fig. 16 shows the 
Monte Carlo histogram with the width of 1.7 mm. 

The temporal stability of the beam stopping distribution, namely its position 
inside the target and its spread, are monitored continuously with the back- 
tracking tomography algorithm. The lateral and the longitudinal centroids 
of the 7r+ stopping spot varied with rms widths of ~ 0.05 mm and 0.2 mm, 
respectively, during the three month calibration period. 

The distances from the forward beam counter BC and the pion production E 
target to the center of the stopping target are 3.87 m and 16.83 m, respectively. 
The e+ and /x"*" times relative to the 7r+ arrival time, calculated under the as- 
sumption of no decay-in-flight contamination, are then —7.3 ns and —2.8 ns for 
the BC-AT path and -12.6 ns and -12.1ns for the E-AT path. Fig. 17 shows 
the beam composition revealed by this method. The expected e'^ / fi'^ TOF 
values quoted above are indicated by the cross markers and agree reasonably 
well with the measured relative timings. The mismatch is due to timing walk 
in the beam counter discriminator. The extracted e"^ and /i'^ beam contam- 
inations in the 7r-in-beam trigger measured in TOF spectra are small, 0.4 % 
and 0.2 %, respectively. 

Calibration runs were performed using a relaxed trigger configuration to study 
the pion beam contamination. We retain a remarkably clean 7r+ beam spot 
with less than 1.0% nonpionic contamination even after omitting the rf signal 
in the tt"*" beam coincidence and reducing the degrader discriminator threshold 
by a factor of three, well below the muon and positron energy depositions. In 
addition, the measured e"*" beam fraction value was confirmed by using an 
extended target with a 50 mm long passive front section that stopped tt"'" and 
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components. Monte Carlo simulations with 7r"'"s generated at the forward 
beam counter position predict an e"*" contamination of 0.5% at the target 
position, arising mostly from /i"*" decays in flight. 



4 Cylindrical MWPCs 



The primary need for a tracking detector in the experiment is due to the 
7re2 decay trigger. Individual positrons from — > euu (Michel) decays can be 
distinguished from the 7re2 decay positrons by their lower energy. However, in 
a high-rate measurement, accidentally coincident Michel events will give rise 
to a substantial background under the 7re2 positron energy peak. 

The solid angle resolution of the calorimeter is ~ 0.004 of 47rsr. This is insuf- 
ficient for the suppression of accidental Michel coincidences occurring at high 
pion stopping rates (> 1 • 10^ /s). In addition to the double-track resolution 
requirement, the following constraints were used in choosing the design of the 
tracking detector: 

(1) low mass, in order to minimize the 7's converting into e'^e~ pairs; 

(2) high efficiency - better than 99.9%; 

(3) high rate capability - up to 10^ minimum-ionizing particles (MIP) per 
second; 

(4) stable operation and good radiation hardness; 

(5) cylindrical geometry. 

After considering different design options we concluded that MWPCs provide 
the best solution. Moreover, wire-chamber technology is mature and reliable. 
We selected a design with a pair of cylindrical chambers, each having one 
anode wire plane along the z direction, and two cathode strip planes in a 
stereoscopic geometry. The MWPCs were manufactured at the Joint Institute 
for Nuclear Research, Dubna. Table 1 lists their basic parameters. Ref. [19] 
gives a general description of the design and operation of the "DUBNA" - 
type cyhndrical chambers. Refs. [20,21] provide a detailed description of the 
PIBETA wire-chambers. In this section we elaborate on the angular resolution 
of the chambers measured with cosmic muons and the chamber detection 
efficiencies measured with minimum ionizing particles. 

The angular resolution of the MWPCs is an important parameter of the 
chamber response. Cosmic muons are ideally suited for the calibration and 
resolution studies of the MWPCs. A clean sample of cosmic muon events is 
collected throughout the experiment during weekly beam-off cyclotron main- 
tenance periods. The trigger logic requires the coincidence of the cosmic muon 
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veto counter (sec Sec. 7) and the "low threshold" (LT~ 5 MeV) Csl calorimeter 
signal (see Sec. 10). 

In off-line data analysis, two additional, stringent cuts are imposed: Csl calorime- 
ter energy deposition more than 200 MeV, and exactly two reconstructed hits 
in each MWPC. These constraints effectively remove any extraneous noncos- 
mic background. The cosmic muon track is reconstructed from a pair of hits 
in one MWPC and the intersections of that track with the other chamber is 
calculated. The difference between the calculated and measured intersection 
coordinates represents the spatial resolution of the chamber. Fig. 18 shows the 
azimuthal angle A0 and longitudinal coordinate Az resolutions of MWPCi, 
properly weighted for the nonuniform angular distribution of cosmic muons. 
These measurements fix the chambers' angular rms resolution at < 0.25°. 

The detection efficiency of the inner chamber for minimum ionizing particles 
eMWPCi can be measured using copious — > e^Ue^iJ. positrons emanating 
from the target: 

_ Ar(AToMWPCioMWPC2oPVoCsI) 
^^^""^^ ~ Ar(AToMWPC2oPVoCsI) ' 



where the A'^'s represent the number of Michel events for which all the detectors 
in the parentheses register coincident hits above the discriminator threshold. 
The Csl calorimeter signal is discriminated with the LT level, while the window 
cut on the PV pulse-height spectrum selects the MIP events. An equivalent ex- 
pression, with indices 1 and 2 exchanged, holds for the outer chamber MWPC2. 
The average MWPC detection efficiencies at the ~ 1 MHz tt"*" beam stopping 
rate are > 94 % and > 97 % for the inner and outer chambers, respectively. 
The combined detection efficiency is therefore > 99.8 %. Fig. 19 demonstrates 
that the experimental charged particle detection efficiency is uniform in the 
azimuthal angle for both chambers. 



5 Plastic veto hodoscope 



The plastic veto hodoscope PV is located in the interior of the calorimeter 
surrounding the two concentric wire chambers. The hodoscope array consists of 
20 independent plastic scintillator staves arranged to form a complete cylinder 
598 mm long with a 132 mm radius. The long axis coincides with the beam 
line and the target axis (see Fig. 20 (i)). The PV hodoscope covers the entire 
geometrical solid angle subtended by the Csl calorimeter as seen from the 
target center. A single PV detector consists of four main components: 
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(1) plastic scintillator stave; 

(2) two light guides; 

(3) two attached photomultiplier tubes; 

(4) aluminized Mylar wrapping. 

The BC-400 scintillator staves were procured from the Bicron Corp. (now 
part of Saint-Goban/Norton Industrial Ceramics Corp.). The maximum light 
emission of the scintillator occurs at a wavelength of 423 nm, while the light 
output falls to 20% of maximum at 410 nm and 480 nm. The rise and decay 
times of the fast scintillator pulses are 0.9 ns and 2.4 ns, respectively. The bulk 
medium attenuation length is 160 cm. The hght response of the scintillator 
has neghgible temperature dependence from — 60°C to +20°C [22]. 

The dimensions of the individual plastic staves are 3.175 x 41.895 x 598 mm^ 
(see Fig. 20 (ii)). They are rigid enough to maintain their shape with support 
only at the two ends. The design geometry ensures that each charged par- 
ticle path emanating from the target volume towards the calorimeter sphere 
intersects the PV array. An 18° angle on one side of each stave minimizes the 
passive PV volumes along the hkely particle trajectories (see Fig. 20(iii)). 

The hght guides flare out to allow more space between the PV PMTs and to 
provide better access to the central part of the detector, facilitating servic- 
ing operations. Each light guide is made of two acrylic pieces that were first 
glued together and then firmly attached to the ends of the PV stave using 
Bicron two-component epoxy glue transparent to the scintillator light. Care 
was taken to eliminate air bubbles in the glue material. After an overnight 
drying period, the glued joints were polished to maintain a precise 3.175 mm 
thickness of the stave. The joint areas were polished to a smooth refiective 
surface by first using a fine grain sandpaper with ~ 600 grit mixed with water 
and subsequently applying a buffing wheel and a polishing gel equivalent to 
10,000-grit sandpaper. The scintillator crazing that can result from the pres- 
sure caused by the weight of the light guides was alleviated by annealing all 
staves in a warm water bath for a period of several hours before installation. 
Each plastic scintillator along with its attached light guides is wrapped in a 
0.25 yum thick aluminized Mylar foil to separate it optically from the adjacent 
staves and to provide the reflective surface which maximizes the amount of 
light reflected towards the phototubes. 

The assembled hodoscope array is supported by the outside surface of a 
530 mm long carbon fiber cylinder with a total thickness of 1 mm, equivalent 
to 5.3 ■ 10~^ radiation lengths. The hodoscope modules are held tightly around 
the support cylinder by a helically wound, thin plastic string, tensioned at the 
extreme ends of the detector stand. 



14 



Before the veto hodoscope was assembled inside the PIBETA calorimeter, the 
light response of each plastic counter was calibrated in a specially designed 
cosmic muon tomography apparatus [23]. We found the average attenuation 
length of the scintillator light to be 396±13mm. The number of photoelec- 
trons per MeV of deposited energy varies between 21 photoelectrons/MeV 
and 63 photoelectrons/MeV, with the average being 38.3 photoelectrons/MeV. 
Fig. 21 shows the measured energy spectrum of the positrons and protons in 
the PV detector, corrected for the angle of incidence. The energy resolution 
measured for minimum ionizing particles is (7^/^=33.2 %. 

The scintillator light is viewed at both detector ends by two Burle Industries 
S83062E photomultiplier tubes. These tubes are 28 mm head-on fast PMTs 
with 10 dynode stages. The photocathodes are sensitive to wavelengths be- 
tween 300 nm and 660 nm. The phototubes are operated at an average voltage 
of 1400 V. 

The individual PMT analog signals from the PV counters are split into two 
branches: one is connected to the trigger electronics and the other is fed into a 
l-to-2 splitter whose appropriately delayed output is wired to the FASTBUS 
ADCs and TDCs. The signals at the trigger side are summed using standard 
NIM logic fan-in modules. The discriminated PV OR signal is used routinely 
as a trigger in MWPC tests. 

The PV array is an essential part of the charged particle tracking system. It 
supplements the MWPC tracking and Csl calorimetry by providing: 

(1) efficient charged particle detection, particularly when combined with the 



(2) reliable discrimination between minimum ionizing particles (cosmic muons, 
positrons/electrons) and protons; 

(3) crude measurements of charged particle azimuthal angle (±9°); 

(4) precise charged particle timing information (±0.8 ns). 

Particle identification is accomplished by cuts on the two-dimensional Epy- 
{E-pY+Ecsi) histogram. This is especially useful for discerning between protons 
and minimum ionizing positrons. Fig. 22 shows the measured two-dimensional 
energy spectra of positrons and protons in the PV detector and Csl calorime- 
ter, corrected for the angle of incidence using MWPC data. The calibrated 
energy boundaries expressed in MeV units are conveniently parameterized by 
the functions 



MWPC data; 



E^y > 0.20 • exp [-0.007 • {Epy + Ecsi)] < E^ 
E% > 2.30 • exp [-0.007 • {Epy + Ecsi)]- 



(2) 
(3) 
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Single positron events fall in a band between the two curves while single proton 
events lie above the second curve. The e+/p identification efficiency Cmp was 
studied in low intensity runs and was determined to be 99.8 ± 0.2 %. 

The PV charged particle detection efficiency is cvahiated separately for mini- 
mum-ionizing positrons with total energies above 5 MeV and nonrelativistic 
protons with kinetic energies in the range 10-150 MeV. A charged particle 
track is defined by coincident hits in the active target AT, the wire cham- 
bers MWPCi and MWPC2, and the Csl calorimeter Csl. The PV detection 
efficiency epv is defined as the following ratio: 

_ Ar(AToMWPCioMWPC2oPVoCsI) 

~ Ar(AToMWPCioMWPC2oCsI) ' ^ ^ 



where the A^'s represent the number of events satisfying the condition in the 
parentheses. The average positron detection efficiency epv measured during 
the detector commissioning period was >99.2% (Figs. 23 and 24). At the 
nominal 1 MHz vr^ stopping rate, the hit rate in an individual plastic scintil- 
lator phototube is ~ 130 kHz while the hit rate in the whole PV hodoscope 
system is ~ 0.88 MHz. Under these conditions, the charged particle tracking 
system, combining the information from the PV hodoscope and MWPC pair 
has a MIP tracking inefficiency in the range (1.0 ± 0.2) • 10~^. 



6 Modular pure Csl calorimeter 



6.1 Calorimeter geometry 



The heart of the PIBETA detector is the shower calorimeter. Its active volume 
is made of pure Cesium Iodide [24-26] . The optical and nuclear properties of 
pure Csl are summarized in Table 2. The experimental signature of a 7r+ — > 
7r°e+i/ event is the prompt decay 7r° — > 77. The calorimeter must be able to 
handle high event rates and cover a large solid angle with high efficiency for 
tt" — s> 77 detection. Efficient suppression of background events requires good 
energy and time resolution. Furthermore, the system must operate with low 
systematic errors and be subject to accurate calibration. The geometry of the 
shower calorimeter is therefore chosen subject to the following conditions: 

(1) high rate-handling ability — low pile-up and low dead time; 

(2) good energy and time resolution for background suppression; 

(3) optimum light collection uniformity for individual modules; 
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(4) uniform solid angle coverage by each module for uniform distribution of 

counting rates; 

(5) maximum economy of different module shapes compatible with the above 
requirements. 

These requirements are best met in the nearly-spherical geometry obtained 
by the ten-frequency class 11 geodesic triangulation of an icosahedron [27]. In 
this respect, our design bears similarities to the SLAC Crystal Ball [28]. How- 
ever, in order to enhance the light collection efficiency, we opted for truncated 
pentagonal, hexagonal, and trapezial pyramids (see Fig. 25) rather than the 
triangular pyramids of the Crystal Ball. 

The PIBETA calorimeter consists of 240 pure Csl crystals (Fig. 26). Our 
geodesic division results in 220 truncated hexagonal and pentagonal pyramids 
covering the total sohd angle of 0.77x47r sr. An additional 20 crystals cover 
two detector openings for the beam entry and exit and act as electromag- 
netic shower leakage vetoes. The inner radius of the calorimeter is 26 cm, and 
the module axial length is 22 cm, corresponding to 12 Csl radiation lengths 
(Xo=1.85cm). There are nine different module shapes: four irregular hexago- 
nal truncated pyramids (we label them HEX-A, HEX-B, HEX-C, and HEX- 
D), one regular pentagon (PENT) and two irregular half- hexagonal truncated 
pyramids (HEX-Dl and HEX-D2), plus two trapezohedrons which function 
as calorimeter vetoes (VET-1 and VET-2). The volumes of the Csl modules 
vary from 797 cm^ (HEX-Dl/2) to 1718 cm^ (HEX-C). Fig. 25 shows aU the 
crystal shapes in nine successive panels. 

The first 25 crystals acquired were manufactured at the Bicron Corporation 
facility in Newbury. Ohio. The remaining 215 were grown in the Institute for 
Single Crystals in Harkov (AMCRYS), Ukraine. Preliminary quality control, 
including the measurement of optical and mechanical properties, was done at 
the factory sites. 

The crystal sphere was designed to be self-supporting. Thus, tight mechanical 
tolerances were required to avoid construction difficulties. Manufacturing tol- 
erances were specified for the linear dimensions (-1-150 /^m/— 50 /xm) and the 
angular deviations (-1-0.040°/— 0.013°). The machined crystals were measured 
upon delivery at PSI using a computer-controlled distance-measuring device 
(WENZEL Precision). The machine was programmed to probe automatically 
the surfaces of the crystal to be tested with a predefined shape. Each crystal 
surface plane was scanned with a touch head at six points which were then 
mathematically fitted with a plane. Body vertices, measured with an absolute 
precision of 20 /im and reproducible to within 2 /im, were then compared with 
the expected values. Crystals that failed the stringent geometrical tolerances 
were returned to the manufacturer for reuse as raw crystal growing material. 
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6.2 Csl detector surface treatment 



After completing the measurements, all crystal surfaces were hand-polished 
with a mixture of 0.2 /xm aluminum oxide powder and ethylenglycol. The 
surfaces of each Csl crystal were then painted with a special organo-silicon 
mixture developed by the Harkov Single Crystals Research Institute [29] . The 
lacquer was a wavelength-shifting ladder organosilicon copolymer with chem- 
ical composition PPO+POPOP+COUM.l, where PPO is 2.5-Diphenyloxa- 
zole, POPOP represents 1.4-Di-2-(5-Phenyloxazolile-Benzene) and COUM.l 
is 7-Diethylamine-4-Methylcoumarin [30] . It provides an optical treatment of 
crystal surfaces superior to more common matting or wrapping treatments. 

The total light output and timing resolution of the painted crystal modules 
were unchanged, but the energy resolution was improved due to the more 
favorable wavelength modulation of axial light collection probabihty. The pro- 
tective surface coating also guards the scintillators against potential surface 
deterioration and minimizes changes of the detector response over the duration 
of the experiment. 

The contribution of the fast decay component to the total light output, the 
photoelectron statistics and the temperature dependence of the light output 
were measured for each Csl crystal using a radioactive source in a specially 
designed scanning apparatus. These measurements arc covered in depth in 
a separate article [31]; here we only summarize the optical properties of Csl 
crystals in Table 2. 

Electron Tubes Inc. (formerly EMI) custom 9822QKB fast 10-stage PMTs [32] 
with 78 mm diameter end windows were glued to the back faces of the hexag- 
onal and pentagonal Csl crystals using a 300 //m layer of silicone Sylgard 
184 elastomer (Dow Corning RTV silicon rubber plus catalyst). The result- 
ing crystal-PMT couplings are strong and permanent, but can be broken by 
application of a substantial tangential force. The smaller half-hexagonal and 
trapezial detector modules were equipped with 46 mm 10-stage EMI 9211QKB 
PMTs [32]. Both photocathodes have quartz windows letting through UV light 
down to 175 nm. The window transparency peaks at ~ 380 nm [32] and is thus 
approximately matched to the spectral excitation of the pure Csl fast scintil- 
lation light component which peaks at ~310 nm at room temperature [33]. 

The PMT high voltage dividers are model UVA 131-a and UVA 131-b, de- 
signed and built at the University of Virginia. The dividers were designed to 

meet several concurrent design goals: (a) high pulse linearity over a signal 
dynamic range exceeding 100 : 1, (b) high gain stability over a wide range of 
counting rates, up to several hundred kHz, and (c) low thermal dissipation 
in order to minimize the heat load on the detector temperature stabilization 



18 



system. In this design the vohagcs of the last four dynode stages were reg- 
ulated using 600 V MOSFET transistors, while the cathode-to-first-dynode 
voltage was fixed by the use of Zener diodes. In this way the dividers reduce 
the so-called "super-linearity" exhibited by many PMTs well below the onset 
of saturation (> 50 jiA average anode current). The choice of SbCs instead of 
the standard BeCu dynode material in the custom 9822QKB tubes further re- 
duces the PMT gain shift that occurs at small anode currents (1-10 /iA). The 
maximum PMT nonlinearity measured in a test with a pair of light-emitting 
diodes was less than 2 % over the full dynamic range encountered in the 7r/3 
decay rate measurement [34] . More details concerning the design, construction 
and performance of the UVA voltage dividers in conjunction with the Electron 
Tubes PMT's will be given in a forthcoming publication. 

Two LeCroy 1440 high voltage mainframes provide the high voltage for the 
PMTs. This 340-channel HV system has computer control capability and thus 
allows for frequent, remote changes in the HV supplied to the PMTs. The 
demand HV values can be set in 1 V increments with reproducibility of ~ 1 V, 
which corresponds to a gain change of < 0.5% for our 10 stage PMTs operating 
in the range of 1500 V to 2200 V. 



6.3 Assembly of the Csl calorimeter sphere 



The Csl calorimeter is mounted inside a forged spherical steel shell which 
is itself supported by a steel frame (the calorimeter mass is approximately 
1,500kg). The shell has two large axial openings, 550mm in diameter, for 
beam entrance and exit, and 220 smaller holes distributed over its surface. 
These smaller holes are aligned with the axes of the individual crystals and are 
intended to allow access to the PMTs (the 20 veto crystals are accessed from a 
neighboring crystal's hole). Inside the shell beam openings are placed cast-iron 
plugs (one on each end) in the shape of pentagonal truncated pyramids. The 
plugs press the stacked crystals together and hold them in place. The plugs 
have a cylindrical borehole (diameter 270 mm) for the beam which also allows 
access to the interior of the Csl sphere, e.g., the wire chambers, the target 
counters, etc.). 

Assembly of the calorimeter proceeded as follows: 

The empty spherical shell was equipped with a single plug insert. The beam 
hole of this plug was used to hold a micrometer gauge with a touch ball 
in place. This device allowed a precise measurement of the distance from a 
crystal's front surface to the center of the sphere. The whole frame was rotated 
by 90° such that the plug was at the bottom. Through the other beam hole. 
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now at the top of the sphere and not yet closed with its plug, all crystals, 
already glued to their PMTs, were moved into the interior of the spherical 
shell. First came 10 pieces of VET-1 and VET-2 type crystals which were 
positioned on the fiat surfaces of the bottom plug. The correct distance to the 
center of the sphere was adjusted using threaded bolts mounted on the plug, 
supporting each crystal's rear face. 

The next crystal layer was stacked on top of the veto layer. Threaded tubes 
with boreholes large enough to house a PMT and //-metal shielding were 
turned through the smaller holes in the shell. The front of these tubes pressed 
against the rear surfaces of the crystals and thereby allowed for a precise 
adjustment of their positions along a radius of the sphere. The high voltage 
divider base was placed in the tube and connected to the PMT. Bases for the 
veto crystals had to be mounted somewhat differently as there are no holes 
available for them on the shell due to the close proximity of the axial plugs. 
The closest hole, later employed for a crystal on the next layer, was used to 
move each veto crystal's base into the sphere, then sideways until connection 
to the veto's PMT was possible. Small supporting structures were mounted 
to lightly press the base against the PMT to prevent loosening due to gravity 
for vertically mounted crystals. 

When the sphere was half filled with crystals, the procedure had to be changed 
somewhat. In order to prevent the crystals in the upper hemisphere from 
sliding down towards the sphere's center, the measurement device at the center 
was replaced by a pin-holder with a pin for each remaining crystal, adjusted 
in length to give the correct distance to the center. 

After the last layer of 10 veto crystals was complete, the second pentagonally 
shaped plug was mounted on top and pressed against the crystal ball under- 
neath. Calculations made prior to the assembly indicated that the plug would 
have to be pushed 10 mm against the crystal ball in order to pack all crystals 
tightly. This prediction was fulfilled exactly, giving us confidence in the correct 
positioning of all modules. The completed sphere was rotated several times 
while being closely monitored for stability before the pin-holder at the center 
was disassembled through the beam entry/exit holes. 



6.4 Calorimeter clustering and nomenclature 



The one-arm calorimeter energy trigger is a basic element of the trigger logic. 
A preliminary simulation study [35] of the calorimeter response to photons 
from ttP decay at rest and 70 MeV positrons from 7re2 decays indicated that: 

(1) electromagnetic shower profiles of the mean deposited energies are similar 
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for photons and positrons, in particular for 6c < 12°, with 6c being a half- 
angle of a conical bin concentric with the direction of an incident particle; 

(2) the average deposited energy and the corresponding energy resolution of 
the calorimeter both reach saturation within a cone of 12° half-angle; 

(3) a centrally hit Csl module receives on average 90 % of the deposited en- 
ergy; at most three modules contain a significant part of shower energy; and 
a group of 9 detectors (a Csl "cluster" ) constitutes an excellent summed en- 
ergy trigger as it registers on average > 98 % of the incident particle energy. 

Therefore, the building blocks of our physics triggers are overlapping clusters of 
nine Csl detectors. Excluding the Csl shower vetoes from the scheme we define 
60 such detector groupings (see Table 3). Every Csl cluster has a symmetric 
partner in the antipodal calorimeter hemisphere. In addition, each Csl module 
belongs to no more than three clusters. This limitation helps to minimize 
the degradation of analog pulses due to excessive signal splittings. In trigger 
design studies that looked at the energy captured by a single cluster as a 
figure of merit, it was found that this clustering scheme, in conjunction with a 
50 MeV discrimination threshold, gives 99.3 % and 98.6 % triggering efficiency 
for 70 MeV photons and positrons, respectively. 

The overall calorimeter topography and the Csl module labels are shown in 
a Mercator projection in Fig. 27. One type of Csl cluster is centered around 
a PENT module (10 pentagons total), and five different types of clusters are 
centered around each of the HEX-B modules (50 hexagon Bs total). Six adja- 
cent Csl clusters are in turn grouped into a Csl "supercluster" : there are ten 
such supcrclusters in the calorimeter each containing 24 individual Csl mod- 
ules. In trigger parlance, a supercluster fires if at least one of its constituent 
clusters fires. A cluster fires if the summed energy of its modules is greater 
than the preset discriminator threshold. 

Ten supercluster labels and crystal types are the basis of the Csl modules' 
nomenclature. The last digit of the crystal number designates the supercluster 
in which a crystal is located. We start by labeling all of them according to their 
central pentagon with numbers 0-9. By separating the calorimeter crystals into 
two hemispheres the beam-upstream crystal labels end in digits 0-4 while the 
names of downstream crystals terminate with digits 5-9. Antipodal crystals 
carry labels differing by 5 units. The crystal map shown in Fig. 27 contains 
module labels as well as the manufacturer's serial numbers (S) of all crystals. 
In this map, the module shape types (P, A, B, C, D, H, and V) precede the 
crystal numbers. 

The serial numbers assigned to the Csl modules run from SOOl to S240: the 
pentagons carry the serial numbers SOOl-SOlO, the hexagon As S011-S060, 
etc. The Csl vetoes receive the last twenty designations S221-S240. The labels 
associated with the clusters and superclusters relate to the actual positions of 
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the modules in the calorimeter, while the serial numbers link the module to the 
online database of crystal properties, such as the shape, machining accuracy, 
fast-to-total light output ratio, photoelectron statistics, light collection unifor- 
mity etc. The same calorimeter database is also used in the GEANT simulation 
of the complete detector response, described below in Sec. 17. 



6.5 Calorimeter electronics 



Each calorimeter PMT is shielded by a metal cylinder that covers both the 
phototube housing and the voltage divider attachment. The voltage divider 
has two anode signal outputs and a high voltage input connector. One analog 
output is delayed in a 75 m (385 ns) long Russian-made PK-50-2-16 coaxial 
cable. These delay cables are smaller in diameter than the more commonly 
used (and more expensive) RG-58 cables (namely 3.5 mm vs 5.0 mm). However, 
they have nearly identical signal propagation speeds, rise times, and peak 
attenuations when tested with a Csl pulse. For a 385 ns delay, they also have 
~ 60 % smaller signal attenuation when compared with RG-58s at 50 MHz and 
100 MHz. 

The delay cables are bundled in groups of sixteen, spooled on six vertical posts 
mounted on the detector platform and connected to the digitizing electronics 
branch. The points of connection to the PIBETA electronics are UVA 139 
ADC signal sphtters that were custom-made at the University of Virginia. 
Each unit is a 16-channel 50 Ohm balanced passive splitter with chip resistors. 
All 16 signal inputs accepting LEMO connectors are provided on the back side 
of the unit, while the front splitter plate features two sets of outputs: a column 
of individual LEMO outputs and a single-header 2x 17 pin connector to couple 
to a twisted ribbon cable. Five UVA 139 units carrying 80 split channels each 
are fitted into a single housing strip with a standard 483 mm electronics rack 
width. 

The LEMO outputs of UVA 139 are routed to CAMAC PhiUips PS 7106 
discriminators with two sets of outputs. One set of discriminated signals is 
taken to a CAMAC Phillips PS 7132H scaler while the other 16-signal twisted 
ribbon output is connected to a FASTBUS LeCroy Model 1877 multi-hit TDC 
unit. The 2x17 simple- header outputs of UVA 139 are split again in the ratio 
1:3 by a custom made plug-in splitter. The split branch with a third of the 
charge is connected to FASTBUS LeCroy Model 1882F ADC, while the parallel 
line is routed to the domino-sampling chip digitizing system (Sec. 9). 

The second set of PMT outputs (identical to the first) are connected to the fast 
trigger electronics section via standard 64 ns long RG-58 coaxial cables which 
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couple to a feedthrough plate where shorter (0.5-10 ns) cables can be inserted 
to equalize the time delays of all Csl signals relative to the time of the slowest 
one. The simultaneous signals are then split in the custom-made UVA 126 3- 
way resistive splitter. This unit has LEMO-type inputs and three outputs per 
input. Each splitter unit contains 60 channels and fits into a standard 483 mm 
wide electronics rack. Up to three copies of Csl analog signals available at the 
outputs arc distributed into dedicated linear summer/discriminator modules 
according to the cluster logic detailed in the previous section. 

The custom-made single- width NIM UVA 125 summer /discriminator is the 
most important trigger unit, performing fast and efficient calorimeter energy 
STimming and discrimination. A single unit has 4 identical sections, each with 
a 9-input summing amplifier followed by a bi-level discriminator. The am- 
plifier output is fed simultaneously to two voltage comparators after passing 
through a protective diode voltage limiting circuit. Each comparator has its 
own reference (threshold) voltage, one low and the other high. The best tim- 
ing is achieved with the low threshold comparator where time slewing versus 
pulse height is less pronounced than it is for the upper threshold comparator. 
The upper level comparator triggers a monostable multivibrator which then 
arms a second monostable multivibrator for a period during which the latter 
can be triggered by the delayed signal from the low level comparator. Thus, 
the low level (LT) and high level (HT) discriminator outputs are essentially 
synchronous. 



7 Cosmic lead house and active cosmic vetoes 



The PIBETA detector is shielded from background radiation in the experi- 
mental hall as well as from the cosmic ray background. The active parts of 
the detector are doubly protected. A lead house enclosure provides the inner 
passive shielding. It is in turn lined on the outside with active cosmic veto 
counters. 

The lead house structure is supported by a steel skeleton frame mounted on 
the detector platform. The skeleton consists of steel beams with 100 x 100 mm^ 
square profiles separated by a 400 mm gap. The upstream side of the support 
structure is closed by a separate lead brick wall fixed to the area fioor. This 
wall has a central hole for the beam vacuum pipe. Successive lead brick rows 
are shifted by half of a brick length starting from the rail pedestal plane. 
Every ~ 300 mm the bricks are attached to the support frame. The lead house 
platform is coupled to the detector platform. Both have a common rail system 
on which the cosmic house can be driven slowly, powered by an electric motor. 
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until it touches the front wall and completes the shielding box which is then 
open only at floor level. A 3 mm thick steel plate is placed on the top of the 
whole structure. 

The active cosmic muon veto consists of 5 extensive scintillator planes shield- 
ing all four sides and the top of the lead house. Each of these planes is made of 
plastic scintillator sandwiched between two green-emitting wavelength shift- 
ing bars, air-gap coupled to its sides (see Fig. 28). The plastic scintillators are 
assembled from 300 mm wide, 25 mm thick panels with lengths of up to 2.29 m. 
They were manufactured at the Scientific Research Concern of the Institute 
for Single Crystals, Harkov, Ukraine. A few scintillator pieces are made smaller 
in order to fit modular planes within the full detector width. The wavelength 
shifter material is BICRON BC-482A extruded into rectangular rods with a 
cross section of 19 x 29 mm^. The cosmic veto counters fining the top and side 
walls are extended 10 cm beyond the walls of the lead house. This geometry 
compensates for the passive volume of the frame structure. All atmospheric 
cosmic muons which hit the PIBETA calorimeter must also intersect at least 
one of the active cosmic veto detectors. 

The optical properties of the plastic counters were tested beforehand and 
found to be comparable with those of BICRON BC-400 plastic scintillator. 

The response of the cosmic veto counter to cosmic muons is illustrated in the 
ADC vs. TDC scatter-plot in Fig. 29. The linearized TDC values are plotted 
along the horizontal axis. The t = point corresponds to the Csl calorimeter 
hit. The ADC readings in noncalibrated channel numbers are shown on the 
vertical axis. Due to the large area of the veto planes and the opportunity for 
the scintillator light to reflect many times in the thin and extensive counters, 
the time resolution is limited to 15 ns FWHM. No timing correction is made to 
account for the varying position of the cosmic muon intersection with the veto 
planes. The energy resolution is consistent with statistics of ~ 5 photoelectrons 
per MeV. The mean detection efficiency is defined as the fraction of cosmic 
events with reconstructed tracks in the calorimeter and the MWPC-PV system 
which also register a TDC hit within a ±15 ns window in at least one cosmic 
veto plane. The mean detection efficiency was measured to be >95 %. 

A cosmic event detected in our apparatus often produces energetic back-to- 
back showers in the calorimeter. If such a cosmic event is in random coincidence 
with a TT-stop trigger and the charged particle detection system is inefficient, 
it might be misinterpreted as a 7r/3 decay event. Fig. 30 (top panel) shows 
the calorimeter energy spectra of cosmic muons taken with a cosmic trigger. 
The trigger is defined as a single Csl cluster firing above the low discrimi- 
nator threshold LT~5 MeV. While the most probable energy deposition is 
~ 170 MeV, calorimeter cosmic energies span the range from the LT thresh- 
old for tracks that just graze the Csl calorimeter, all the way up to almost 
1000 MeV for the extended cosmic showers. 
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The calorimeter 7r/5 trigger requires two coincident superclusters to be located 
in the opposite calorimeter hemispheres firing above the high discriminator 
threshold HT~52MeV. The measured sum of the two cluster energies for 
16, 916 cosmic triggers (a single 1.3 hour long cosmic run) shown in the bottom 
panel of Fig. 30 has a low energy cutoff at ~ 200 MeV. We can use this value as 
an upper level cut in the offline nf] event analysis and thus achieve a minimum 
suppression factor of 2.3 ■ 10^^. The raw counting rate of HT cosmic triggers is 
~ 5 Hz while the nominal tt"*" stopping rate < 1 MHz. With logic signals about 
10 ns wide, the random cosmic-TT-stop overlap rate is 5-10~^Hz. 

Combining the calorimeter energy cut defined above and the 1 ■ 10^^ track- 
ing system inefficiency discussed previously we find that the cosmic back- 
ground contributes to the 7r/3 event sample at a level of <10^'^. Additionally, 
information from the cosmic muon veto counters, which enclose the detec- 
tor, suppresses this contamination by another order of magnitude, making it 
negligible. 



8 Temperature stabilization control 



The design goals for the temperature stabilization system were: 

(1) constant Csl calorimeter temperature of 22° C stable to ±0.2° C; 

(2) relative humidity inside the detector thermal housing < 50%. 

If the air is cooled in the heat exchanger to just above the freezing point, 
the humidity condenses and the water is removed from the thermal volume 
through a drain pipe. A temperature of 0.2° C gives a vapor pressure of 460 Pa, 
corresponding to a relative humidity of 25% at 22°. With this method, dry 
air is produced, and the water sensitive Csl crystal surfaces are protected. 

We can calculate the cooling power requirement as follows: the surface area 
of the thermal enclosure is 28 m^. The Styrodur 4000S panels used for the 
thermal enclosure have a heat conductivity of 0.028 W/m-K. Assuming an 
extreme operating temperature of 10° C with an outside temperature of 35° C 
and insulator thickness of 0.04 m we find a maximum heat load of 

_ 28m2-25K- 0.028 W 

P = = 490W. 5 

0.04m-m-K ^ ^ 



This power dissipation has to be increased by another 500-600 W to account 
for the detector PMT bases, and 60 W for the air circulation blowers. A HX- 
200 Recirculating Chiller from NESLAB Instruments [36] suffices and allows 



25 



some margin of error. It is able to dissipate ~2000 W at 0°C. This unit is 
designed to maintain the temperature of the chilled liquid (usually water) 
to within ±0.1° C. The chiller unit can also be regulated from an external 
temperature sensor. It requires an analog input signal and produces an analog 
output signal. The cooling fluid is kept at a constant temperature of 0.2°C 
using the regulation of the HX-200 unit while a second feedback loop is used 
to regulate the temperature of the air in the enclosure with a heater in the air 
stream. 

Since the temperature is above its freezing point, water (with some antifreeze 
for safety) is an acceptable cooling fluid. The high-flow centrifugal pump (such 
as the PD2 Model from Sta-Rite Corp. [37]) can handle 11.3 i/min which pro- 
vides a temperature differential of 1.5° C. To avoid potential erosion problems, 
two heat exchangers in parallel, both running at about 10 i/min are used giv- 
ing a temperature difference of 0.85°C in the coohng water. This configuration 
also has the advantage of permitting cold air to be blown in from two sides. 

Two heaters, each with a heating power of 2.1 kW, are used in the air pipes 
just after the cooling heat exchanger to keep the temperature stable under 
varying conditions. A regulator turns the heaters on and off using a Triac 
circuit with a period of about one second. Switching happens at the zero 
crossing of the line power, in order to reduce the influence on the trigger 
and data acquisition electronics. The duty cycle is determined by an external 
voltage in the range 0-5 V, which can be either generated manually with a 
potentiometer or controlled by a computer. 

The air temperature is measured just after both heaters, and the Csl temper- 
atures are measured with eight sensors glued directly on the front and back 
faces of Csl crystals at different positions. The slow control computer imple- 
ments two proportional integral-differential (PID) loops in software. The outer 
loop tries to keep the Csl temperature constant by means of varying the tem- 
perature of the heated air. The time constant of this loop was experimentally 
determined to be about one hour due to the large mass of the Csl crystals and 
the metal frame. The inner loop receives the demand value for the air tem- 
perature from the outer loop, and controls the heating power in order to keep 
the air temperature at the demand value. The time constant for this loop is 
determined by the mass of the heating coils and was set at about 30 seconds. 

The combination of the two PID loops keeps the temperature of the Csl crys- 
tals extremely stable. The stability is demonstrated in a four-day snapshot 
of the Csl temperature, averaged over 6 sensors (Fig. 31). The sensor values 
are sampled several times each second and sent to the history subsystem of 
the data acquisition system whenever they change significantly. There they 
are written to a hard disk. The temperature oscillations have an rms value of 
0.02° C. They are due to day-night changes of the outside temperature and 
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do not cause detectable Csl light output variations. During the same time pe- 
riod the relative humidity inside the detector enclosure was regulated within 
a ~ ±3 % range. 



9 500 MHz Domino Sampling Chip 



The Domino Sampling Chip (DSC) is a significant addition to our electronics 
arsenal. It is used to digitize waveforms from every phototube in the detec- 
tor [38]. Waveform digitization helps achieve higher background suppression 
and thus, improved calorimeter energy resolution at higher event rates. Just 
as important is the improvement of time resolution. Our current analysis in- 
dicates that 200 ps rms accuracy in the relative time of the leading edge of 
digitized Csl pulses is attainable. The ultimate accuracy of the method de- 
pends on the noise level in the actual signal pulse shapes. 

The DSC was originally considered for a more advanced stage of the experi- 
ment. However, successful initial tests made the chip available sooner. A ded- 
icated PIBETA circuit board with zero suppression and appropriate readout 
features was developed at PSI. The device is now fully integrated into the PI- 
BETA apparatus and routinely operational. In addition to its other benefits, 
the DSC is also fast (up to 1.2 GHz) and inexpensive. 



9.1 DSC scale calibration 



The DSC was first used with success, albeit with only 5 beamlinc detectors 
(BC, ACs, AD, and AT), during detector commissioning in 1999. Full im- 
plementation, with digitization of all detector channels, was attained in the 
summer of 2000. Each chip operates at ~500MHz producing waveforms 256 ns 
wide in steps of about 2 ns. The analog waveforms are pre-analyzed in the DAQ 
front-end computer. Zero suppression is applied in software to reduce the data 
rate by a factor of ten, while retaining the interesting sections of the waveforms 
for offline analysis. 

The domino speed and thus the sampling frequency varies from channel to 
channel. Though all channels have been set to ^ 500 MHz by manually ad- 
justing the voltage levels which control the domino speed, a precise calibra- 
tion is accomplished in software, taking care of the differences in the preset 
voltages as well as temporal time scale drifts. The random trigger induces a 
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100 ns calibration signal in order to determine the exact sampling speed. The 
signal, capacitively coupled to each DSC input, consists of two voltage spikes 
set exactly 100ns apart. The exact positions of the spikes are determined by 
fitting a second-order polynomial through three peak points, resulting in sub- 
bin time resolution. DSC speed in a specific run is determined by analyzing 
and averaging each channel over all random triggers. Since the random trig- 
ger events come at ~ 1 Hz frequency interlaced with our physics triggers, the 
calibration process runs permanently in the analyzer and updates the DSC 
calibration constants about once an hour. The examination of DSC calibra- 
tion runs taken one week apart testifies to the stability of the domino speed. 
The maximum run-to-run difference is 0.3% which corresponds to 0.4 ns for 
an average peak in the center of the time scale. 



9.2 DSC waveform fitting 



A digitized PMT waveform can be decomposed into a collection of ("ADC",- 
"TDC") pairs corresponding to the individual signals present. The waveform 
analysis thus reduces to determination of the exact shape of an average wave- 
form (a "system function" g) for each particle type (pion, positron, photon or 
proton) in each PMT channel, followed by fitting of all waveforms with these 
functions. If there is more than one hit in a given 256 ns wide window, the 
analyzer uses a sum of system functions, each with a different multiplier 
(related to the ADC values) and adjustable time offsets (corresponding to the 
TDC values): 

N 

fit) = ^DADQ • g{t - DTDQ), (6) 

i=l 

where g{t) is a system function for a given PMT channel and f{t) is the 
function being fitted to a waveform by varying DADCj and DTDCj values in 
order to minimize the of the fit. 

Fig. 32 shows a typical system function for a vr"*" in the active degrader counter, 
derived by averaging waveforms of ~ 100,000 7r+-beam events which contained 
only one hit in that particular channel. Ripples such as the one around t — 
175 ns are due to either refiections on imperfect 50 Ohm terminations in the 
signal path, or to the delayed arrival of positive ions inside the PMT. In either 
case, their amplitude is proportional to the main signal amplitude. These 
ripples are present in this channel for every event and are correctly taken care 
of by the system function. Other techniques, such as using multi-hit TDCs, 
would falsely identify this peak as a real hit. The system functions for all DSC 
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channels arc stored in a waveform database file which resides in the analyzer 
directory and must be present to run the code successfully. 

The single event display in the Fig. 32 shows a representative 7r+ waveform 
from the active degrader channel (dark line) with the fitted function super- 
imposed (lighter line). The flat waveform sections arise from the front-end 
zero suppression, which stores only the data around real hits. The sections 
with small amplitude and small curvature are not taken into account in fit- 
ting. Open circles indicate the extracted DTDC values. The time axis on the 
presented plot is not calibrated, but one horizontal unit corresponds to ~ 2 ns. 

Comparison can be made between the data derived from the waveform dig- 
itizer and traditional FASTBUS and CAMAC ADC/TDC values since we 
acquire both data streams in parallel for all PMT signals. Four panels of 
Fig. 33 show correlation plots between selected variables. The top left panel 
shows the DTDC values plotted against the FASTBUS TDC values. The gap 
around zero is caused by the trigger veto, which suppresses prompt events. To 
make a meaningful comparison with the ADC data, further cuts are helpful. 
The second panel (top right) shows raw DSC data on the vertical scale vs. 
CAMAC ADC data plotted on the horizontal scale. Since the width of the 
ADC gate is only 25 ns, analog signals that are not in time with the calorime- 
ter trigger are not completely digitized during a CAMAC ADC gate, and 
therefore show smaller ADC values. This limitation is not present in the DSC 
data, which covers the full range of 256 ns. One can see on the same panel 
a distribution which, when projected onto the vertical axis, has a peak near 
channel 270 corresponding to a valid pion signal. The projection onto the hor- 
izontal axis produces only a broad distribution extending all the way to zero 
time, corresponding to events where the incoming pion comes > 25 ns before 
the trigger and is therefore absent in the CAMAC ADC data. Only the band 
along the main diagonal shows a clean correlation between DSC and ADC 
data. Selecting prompt trigger events, as in the third panel, only this band is 
present. The correlation arises from the fact that prompt events signals always 
come at t = ns and are therefore fully contained inside the ADC gate. The 
last panel of Fig. 33 shows DSC ADC vs. DSC TDC data. There is no visible 
dependency, as expected, unlike for the CAMAC ADC data. 

A closer look at the time resolution reveals the 0.5 ns resolution of the LRS 
1877 TDC; there is also a timing jitter between the DSC and TDC variables. 
The differences between DTDCs and FASTBUS TDCs have a distribution 
with 0.4 ns rms which is dominated by the limited LRS 1877 TDC resolution. 
A more detailed presentation of the PIBETA waveform digitizer system will 
be published in a separate paper. 
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10 PIBETA detector fast trigger logic 



Selective, bias-free triggers capable of handling high event rates are an essential 
requirement of the detector system. Fast analog triggers were optimized to 
accept simultaneously: 

(1) beam particle events; 

(2) prompt interaction events; 

(3) delayed pion and muon decay processes; 

(4) cosmic muon background; 

(5) true random calibration events. 

A complete list of the 12 triggers used is given in Table 4. The table summarizes 
trigger names, trigger logic definitions, prescaling factors (if any), and typical 
raw event rates recorded at the production 7r+ beam fiux of < IMHz. Here we 
discuss the predefined trigger classes in order of increasing their complexity. 



10.1 Random (ADC pedestal) trigger 



A 190x20x8 mm^ plastic scintillator counter made of BICRON BC-400 mate- 
rial is placed above the electronics racks, parallel to the area fioor, about 3 m 
away from the main PIBETA detector. By virtue of its position, the counter 
is shielded from the experimental radiation area by a 50 mm thick lead brick 
wall as well as a 500 mm thick concrete wall. Operating with a high discrimi- 
nator threshold, it counts only cosmic muons and random background events 
at about l-2sec~^, and has a stable counting rate independent of the beam. 
Its discriminated signals define our true random trigger. 

The random trigger Tr is connected to the "trigger mixer" logic unit and 
during production runs is always included among the enabled triggers. A sin- 
gle production run is limited to 200,000 events, which includes ~10^ random 
trigger events. The onhne analysis program updates the ADC pedestal spec- 
tra of all ADC channels for every random event. The pedestal peaks are then 
analyzed and realigned automatically by an end-of-run routine. A detailed 
description of the ADC pedestal correction algorithm is given in an earlier 
paper [39]. 

We note parenthetically that in off-line data analysis, the ADC values for all 
detectors in random events are written to separate files. These random event 
energy depositions are used subsequently in a GEANT simulation [11] of the 
PIBETA detector, enabhng us to account rigorously for the residual ADC 



30 



noise and accidental event pile-ups in the energy spectra. Fig. 34 shows the 
calibrated calorimeter energy spectrum for the random trigger events. The 
average deposited energy per single Csl module associated with the random 
events at our production tt"*" stop rate is 0.15 MeV. The random trigger comes 
most often in accidental coincidence with a Michel decay event. Therefore, the 
random calorimeter spectrum in the figure has a steep cut-off corresponding 
to the ~ 53 MeV Michel end-point energy minus ~ 4 MeV energy losses in the 
target and PV counters. 



10.2 Beam particle triggers 



The particle-in-beam trigger Tb is defined by a four-fold coincidence between 
the forward beam counter BC, the active degrader AD, the active target AT, 
and the rf accelerator signal: 

Tb = BCoADoATorf. (7) 

Minimum-ionizing positrons in the BC, AD, and AT counters deposit 0.6 MeV, 
7.2 MeV, and 9.0 MeV, respectively. The corresponding energy depositions for 
114.0 MeV/c pions (with kinetic energy 40.6 MeV) are 0.7 MeV, 12.7 MeV, 
and 28.0 MeV. By appropriately adjusting the discriminator thresholds and 
the relative timing of inputs into the quadruple coincidence (7), the vr-in-beam 
and e-in-beam triggers T^r and Tg are set up. Fig. 6 shows a digital oscilloscope 
snap-shot of the four input signals constituting a 7r-in-beam coincidence. 

Moreover, we define a less restrictive 7r-beam signal B^ that tags a t:^ beam 
particle inside the detector via an overlap between the forward beam counter 
pulse BC and the active degrader pulse AD: 

B^ = BCoAD. (8) 

Discriminator thresholds and relative times in the beam particle coincidence 
Tb are adjusted to preferentially select through-going pions. 

Each pion stopping in the target initiates a pion gate vrG, a 180 ns long window, 
whose delay is adjusted to start 50 ns ahead of the pion stop time to- A 10 ns 
wide gate ttS coincident with a pion stop pulse is also generated. Tb, ttG, and 
ttS are used as the building blocks of more complex triggers described below. 
Higher order triggers are generated on the basis of a coincidence (or anti- 
coincidence) between one of these gates and shower signal(s) in the calorimeter. 
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10.3 One-arm calorimeter trigger 



The 15 UVA 125 discriminator/summers define 60 discriminated Csl cluster 
pulses and provide both low and high threshold discriminated outputs. Fig. 35 
shows a schematic layout of the fast trigger logic beginning at the discriminator 
outputs. Two CAMAC LeCroy Model 4564 OR logic units are used to combine 
these 60 clusters into 10 Csl supercluster pulses. Again, both high (HT) and 
low (LT) supercluster variants of the output signals are implemented. The 
simplest calorimeter trigger is the one-arm Csl trigger Cs. It requires the 
firing of at least one Csl supercluster: 



where the numbers 0-9 represent ten Csl superclusters. The OR logic of Eq. (9) 
is implemented in a single PSI LB 500 logic box, described in the next section. 

7re2 decays, for example, are identified by means of a prescaled trigger that re- 
quires detection of a single localized calorimeter shower with deposited energy 
exceeding the HT level during a prescaled pion gate ttGps: 



where 6,^ represents a 10 ns wide prompt veto signal overlapping the vr stop time. 
This coincidence logic is illustrated in Fig. 36 which shows a digital oscilloscope 
snap-shot of the trigger inputs. 



10.4 Two-arm calorimeter trigger 



The default vr/? event calorimeter logic requires a coincidence between any Csl super- 
cluster and a logic OR of five superclusters that are not its neighbors. The bottom 
panel of Fig. 26 shows the calorimeter hemisphere complementary to a single Csl 
supercluster. This relaxed requirement simultaneously samples a broad unbiased 
non-7r/5 background. The logic requirement could be made more stringent for higher 
running rates. However, the experiment also has the goal of measuring radiative 
pion and muon decays with broader angular correlation coverage and this requires 
the logic as implemented. 

In terms of the Csl supercluster indices, the two-arm calorimeter trigger Cgg is 
defined as: 



Css = 0o(2 + 3 + 5 + 6 + 9) + lo(3 + 4 + 5 + 6 + 7)+ 2o(0 + 4 + 6 + 7 + 8) + 



Cs = 0+1 + 2 + 3 + 4 + 5 + 6 + 7 + 8 + 9, 



(9) 



T^ = C^oTrGpgoB. 



(10) 
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3o(0 + 1 + 7 + 8 + 9) + 4o(l + 2 + 5 + 8 + 9) + 

5o(0 + 1 + 4 + 7 + 8)+ 6o(0 + 1 + 2 + 8 + 9)+ 7o(l + 2 + 3 + 5 + 9) + 
8o(2 + 3 + 4 + 5 + 6) + 9o(0 + 3 + 4 + 6 + 7), (11) 



where the first two hnes refer to the "northern" calorimeter hemisphere, and the 
last two lines refer to the "southern" hemisphere. The trigger OR involves the logic 
sum of 10 conjugate ORs that are implemented in the first programmable LB 500 
logic box. The actual logic function implemented omits the repetition of two-fold 
AND terms. 

The two-arm trigger has three versions, one with high discriminator thresholds on 
both superclusters Cgg, another with two low discriminator thresholds Cgg, and a 
third version with one high and one low threshold Cg|^. The latter versions allow 
sampling of the radiative tt and /x decays. 

7r(3 decays as well as radiative vr — > ei"y events are identified by means of a dedicated 
trigger requiring two antipodal Csl showers during the ttG with energy deposited 
in the calorimeter exceeding HT: 

Tss = CiovrGoB,. (12) 

The radiative muon decays — > euv^ are captured with a prescaled two-arm low 
threshold calorimeter trigger: 

Tss = C^s°^GoB,. (13) 

Finally, the third version of the two-arm trigger samples dominantly the radiative 
TT — > euj decays with a low energy positron or photon: 

Tsf = CLgSoTTGoB,. (14) 



10.5 Three-arm calorimeter trigger 



The three-arm calorimeter trigger C3S is defined as: 

C3S = 0o2o9 + 0o3o7 + 0o7o9 + lo3o5 + lo4o8 + 

lo5o8 + 2o4o6 + 2o6o9 + 3o5o7 + 4o6o8. (15) 

The ten terms in the logic sum represent all possible combinations of three nonad- 
jacent low-threshold supercluster overlaps. Like the others, the three-arm trigger is 
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implemented in the LB 500 unit. We use this trigger to eollect events with three 
coincident particles in the final state following the tt stop signal: 

T3S = CasoTTGoB^. (16) 

This trigger is dominated by the random coincidences of three Michel positrons but 
it also includes the Dalitz decays ir^ — ^ e'^e~7, and rare fi'^ e^v^i'ne^e' decays. 



10.6 Prompt event trigger 



The prompt event trigger Tp uses the one-arm Csl logic Cs and an adequately 
prescaled version of the pion stop TrSps' : 

Tp = Cf ottSps'. (17) 

During production running we select a high threshold version of the Cs input, 
whereas the energy spectrum of prompt protons in the calorimeter shown in Fig. 37 
is taken in a calibration run using the low-threshold version in order to record the 
entire energy spectrum of prompt events. 



10.7 Cosmic muon event trigger 



There are two different versions of the cosmic muon trigger. In routine production 
running with high intensity Tr"*" beam, a cosmic event is defined as a coincidence 
between a prescaled signal from the cosmic muon veto counter CVps and a high 
threshold one-arm Csl signal Cg , vetoed by the beam: 

tH = CVps o C^ o B,. (18) 

The low-threshold one-arm Csl version of the cosmic trigger T^ defined as 

tL = CVoC^oB^, (19) 

uses a nonprescaled cosmic muon counter signal CV. It is included in the trigger 

mix only during periods where there is no tt~^ beam in order to take cosmic muon 
data for calibration purposes. This usually happened during planned accelerator 
maintenance - normally one day each week. 
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10.8 The PSI LB 500 programmable logic unit 



The several triggers described in the previous sections were originaUy implemented 
in traditional CAMAC coincidence logic units with specially made twisted pair ca- 
bles. This solution proved to be unreliable and inflexible. Therefore, a new logic 
unit denoted LB 500 was developed at PSI. This single width CAMAC unit con- 
tains ECL-TTL converters for 64 input/output lines and a Lattice Complex Pro- 
grammable Logic Device (CPLD) ispLSI 2128E [40]. This chip has 6000 programmable 
gates which can be combined in virtually any way and a pin-to-pin propagation de- 
lay of 5 ns. Unlike field programmable gate arrays (FPGAs), the CPLD does not 
contain static RAM cells for logic functions. Therefore, it can function exactly as 
a conventional coincidence unit in overlap mode, preserving the times of the input 
signals at the output. The LB 500 unit can be reprogrammed via CAMAC in a few 
seconds. This allows remote modification of the trigger in a reliable way. The 64 
ECL lines at the front panel of the LB 500 unit can be configured in groups of four 
as either inputs or outputs by changing internal resistors. 

One LB 500 unit is used for defining the two-arm and three-arm triggers and another 
unit is used for prescaling selected triggers. This solution is more powerful than that 
using conventional prescalers with clocks and coincidence units, since the prescaling 
factor can be remotely adjusted over a wide range. Internally, the unit uses 15-bit 
synchronous counters, which have propagation delays independent of the prescaling 
factor. Reshaping of the output signal preserves the original input pulse width and 
time. 

A third LB 500 unit is used for the trigger mix. The signal timing inside the LB 
500 unit can only be adjusted in steps of ~2.5ns by means of inverter chains. An 
external programmable CAMAC delay LeCroy Model 4418 is used to adjust the 
timing of the calorimeter triggers and the various beam gates with an accuracy of 
1 ns. In normal production mode, all predefined triggers are enabled except the low 
threshold cosmic trigger. 

Fig. 35 referenced above is a fiow chart of the signal processing described above. 
The design details for all three LB 500 units are available on the World- Wide- Web 
at www.pibeta.psi or www.pibeta.phys.virginia.edu. 



11 The MIDAS data acquisition system 



A new, general purpose data acquisition system called MIDAS (Maximum Inte- 
grated Data Acquisition System) [41] was developed for the PIBETA experiment. 
It is now used in a number of experiments at several laboratories. The system is 
suitable for experiments ranging from small test set-ups with only one personal com- 
puter (PC) connected to CAMAC, through medium scale experiments with several 
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front-end and backend computers. The system runs under Linux, Windows NT, and 
Vx Works. 

An integrated slow control system includes a fast online database and a history 
system. Drivers exist for CAM AC, VME, FASTBUS, high voltage systems, GPIB 
and several PC plug-in DAQ boards. A framework is supplied which can be extended 
by user code for a front-end readout on the one side, and data analysis on the other. 
The online data can be presented in real time via the CERN Physics Analysis 
Workstation library (PAW) in histogram and N-tuple format [42]. The run control is 
done via a Web interface, which allows remote monitoring through the Internet. 

MIDAS software consists of a set of library functions and applications. The library 
is designed in two layers. The lower layer covers all calls dependent on the operating 
system, mainly shared memory access and semaphore operations. The upper layer 
is entirely independent of the operating system. Since the lower layer contains only 
about 7 % of the code, porting the library to a different operating system is quite 
easy. 

The MIDAS library contains routines for buffer management, a message system, 
a history system and an online database (ODE). MIDAS buffers are FIFOs which 
support multiple producers and multiple consumers. Consumers may request a cer- 
tain subset of event types from a buffer. They can also specify whether they want to 
receive all events of a given type or only as many as they can process without block- 
ing the producers. A watchdog scheme has been designed which removes crashed 
clients from a buffer to prevent them from blocking the whole system. The transfer 
speed between a producer and a consumer is on the order of ~20MB/s if both 
run on the same computer and ~ 10 MB/s over 100BaseT Ethernet. The history 
system is used to store data on a hard disk and produce value-versus-time plots. It 
is capable of storing several thousand events per second and of changing the event 
definition on-the-fly during an experiment. 

The online database provides central data storage and contains all relevant experi- 
ment variables like logging channel information, event definitions, slow control vari- 
ables, front-end parameters and histogram definitions. It is kept entirely in shared 
memory for fast access of up to 50,000 read/writes per second locally and 500 
read/writes per second remotely. The database is hierarchically structured, simi- 
lar to a file system, with directories and sub-directories. The data arc stored in 
key /value pairs, where a key corresponds to a file name and the value to the file 
contents. Keys can reside in directories which can be subdirectories themselves. 
ODB keys can be created, edited and deleted dynamically during run-time. Each 
key can contain either a single value of any type (an integer, real, string, etc.) or an 
array of values of the same type. 

Entries in the ODB can be changed in three ways. The first possibility is direct 
editing using ODBEdit, a general purpose editor developed to view and change 
values in the ODB. This method is used mostly when running online. Parameters 
can be changed and the effect of the change can be inspected immediately by look- 
ing at the online N-tuples and histograms. Access to the online database is also 
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possible through a Web interface. Fig. 38 shows an example of the main PIBETA 
experiment control window. This makes it possible to control an experiment from 
any computer running a Web browser. The final way to change ODB variables is 
by loading configuration files. Subtrees of the ODB can be saved and loaded in a 
simple ASCII format. 

A client can register a "hot-link" of a local C-structure to a database sub-tree. 
Whenever a value in the sub-tree is changed by someone, the chent C-structure 
automatically receives an update of the modified parameter. This scheme makes it 
very easy to control dynamically the behavior of front-end programs and analyzers. 
The slow control system is based on the same principle. Whenever a demand value 
is changed in the online database, the slow control front-end receives an update 
and can propagate it to the proper hardware via a device driver. An image of the 
whole database contents can be written to the logging channels to reflect the current 
status of the experiment. 

Several applications have been written Tising the MIDAS library. A general-purpose 
logger supports multiple logging channels to disks, magnetic tapes and ftp sites. 
Different logging channels can receive different types of events in different formats. 
Three data formats are currently supported: ASCII, an optimized MIDAS format 
and the YBOS [43] format. 

A general-purpose analyzer framework was designed which incorporates CERN 
HBOOK library routines [44]. Raw and derived data values are booked automatically 
as N-tuples and can be histogrammed and presented with PAW commands online. 
The identical analyzer can be used for an offline replay analysis. 

The MIDAS implementation used by the PIBETA experiment has two front-end 
PC computers connected to CAMAC, FASTBUS, VME and a LeCroy LRS 1440 
high voltage system. The event response of the first front-end computer running 
under Windows NT using a polling scheme is 8 //s. The second front-end computer 
measures temperature, humidity and gas flow values of the MWPC system and 
controls the high voltage and detector temperature as described in Sec. 8. 

The backend computer running under Linux performs a full online data analysis 
and stores data on disk and magnetic tape. Data files are written to a dual digital 
linear tape drive (DLT TTi 2200) with 20 GB capacity and 1.5 MB/s transfer rate. 
In addition, the data files are written to a disk, from which they are migrated 
automatically to the PSI file archive. The archive system uses a HP/CONVEX 
SPPIOOO mainframe with several DLT tape robots and UniTREE software. Data 
can be transferred using the FTP protocol with a dedicated Ethernet connection 
from the PIBETA experiment. Logically, a user can think of the backup archive 
simply as another disk directory tree. From a user perspective, archive files appear 
as if they were on disk. All 3.4 TB of data recorded by the PIBETA experiment 
could be replayed directly from the archive. 

The MIDAS source code for Unix (Digital Unix/Ultrix/Linux) and the executa- 
bles for Windows NT (Intel) and MS-DOS are available at the MIDAS http server 
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http : //midas .psi . ch/ or http : //midas . triumf . ca. The Concurrent Version Sys- 
tem (CVS) tree that lets developers access the latest code is available via the cvsweb 
interface. 



12 User DAQ code 



12.1 PIBETA analyzer 

The PIBETA analyzer program is written in the C language. It consists of two 
parts: a system part which is responsible for reading and writing events in various 
formats and a user part which actually does the experiment-specific data analysis. 
The analyzer format is modular and user-friendly. It is easy to install, use and 
extend. 

In order to make the data analysis more flexible, a multi-stage concept is chosen as 
the backbone of the analyzer structure. A raw event is passed through several stages 
in the analyzer, where each stage has a specific task. The different stages read an 
event, analyze it and append the results of the analysis back to the event. Therefore, 
each stage in the chain has access to all results from previous tasks. The first stages 
in a chain typically deal with data calibration, while the last stages contain the 
"physics" analysis code, which calculates kincmatical and physical variables. The 
multi-stage concept allows the collaboration members to use standard modules for 
the calibration stages ensuring that all experimentalists deal with identically cali- 
brated data, while the last stages can be modified by individuals to look at different 
aspects of the data. 

This method is different from the usage of data summary tapes (DSTs) in other 
experiments. Instead of producing some intermediate data which get distributed, 
analysis is always performed on the original raw data. The advantage of this method 
is that one still has access to the rudimentary data (such as raw TDC and ADC 
values) even in the most advanced stages of the analysis. If one is in doubt about the 
calibration, one can always go back and test a different calibration method, which 
is not possible if one works on precalibrated DST data. 

A bank system used for event storage reflects the multi-stage concept in the data 
structures. A bank contains a sub-structure of an event. It can be of variable length 
(like the sparsified ADC data) or of fixed length. A fixed size bank can contain 
different data types like integers and floating point values, while a variable length 
bank can only contain a single data type. The front-end program produces banks 
related to different parts of the detector electronics, like a raw ADC bank, a raw 
TDC bank, a scaler bank and so on. The first analyzer stages use these banks to 
produce calibrated data such as energy deposition expressed in MeV, residing in 
a calibrated ADC bank, a charged particle track intersection with the horizonatal 
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plane in millimeters in a MWPC bank and so on. Different users can add private 
banks wfiich contain variables they are interested in. 

Several analyzer modules use banks from the online system to produce new, calcu- 
lated banks. At the end of an analysis process, all events are written to disk. Each 

bank has a database flag telling the system if this bank should be included in the 
output file. By suppressing some extensive online banks, the amount of data in the 
output file can be reduced significantly. 

Since the contents of banks are defined in the online database, the system part 
of the analyzer knows how to interperet the contents of an event. Thus, N-tuples 
can be booked automatically by the system. When running the analyzer offline, 
column-wise N-tuples (CWNTs) are used. Each event bank is booked as a "block" 
in a CWN-tuple. The CWNT identifier is identical to the event number. 

At present the PIBETA analysis at each analyzer stage is performed by 16 indepen- 
dent modules. Each module is written in a different source file and exports routines 
which are called for each event at the beginning and end of a run. Currently, the 
modules have to be written in the C language, though MIDAS does permit FOR- 
TRAN analyzers. Each module uses a set of parameters, which are stored in the 
online database. The advantage of this scheme is that the parameters can be changed 
without recompilation of the analyzer. 

Modules can be placed into cither standard modules or private (user) modules. 
Standard modules are generally accepted by all collaboration members and used in 
every analysis pass, e.g., modules for energy and time calibration, ensuring uniform 
calibration for all users. User modules are written by individuals to look into specific 
aspects of the data. Over time, user modules which are accepted by the collaboration 
can be promoted to the standard modules. If a specific method is established during 
offline analysis, it can go straight to the online analysis of the next beam time, so 
one gets online the same results that previously came only after offline analysis. 

The online database stores all variables which concern a specific experiment. The 
same database is used both online and offline. When running the analyzer offline, 
user configuration files can be loaded, overwriting the ODB parameters stored in the 
raw data file. Thus, different configuration parameters can be loaded for different 
individual runs or different series of runs. 

The analyzer task can read data from hard-disk files, tapes or the online DAQ 
system. It can write events or specific event banks to the output files. The program 
currently supports three different output file formats: MIDAS binary, pure ASCII 
files, and HBOOK RZ files with both Column- Wise and Row- Wise N-tuples. While 
the RZ file format can only be used for the analyzer output, binary and text files 
are readable by the analyzer. Moreover, these types of files can be written and read 
directly in GNU-zipped format: the data stream is compressed and decompressed 
on-the-fly. While this method reduces disk space requirements by 50 %, it takes 
about 20 % more CPU time. 
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The following online run types are predefined for the PIBETA experiment: 

(1) physics data acquisition run; 

(2) pedestal run; 

(3) timing run; 

(4) cosmic run. 

The pedestal run is performed when bringing the detector online after long oper- 
ational interruptions. Its routines update the ADC pedestal values for all detector 
channels. A run in the timing mode should follow by checking the trigger times and 
alerts the user if any adjustments in the time delays of individual sub-detectors are 
necessary. During production running, the physics DAQ mode is chosen for stable 
beam-on periods while the cosmic running mode is used during longer beam-off or 
low beam intervals. Depending on the run mode, the analyzer books different his- 
tograms optimised for that run mode. Since the analyzer has access to the whole 
hardware through the ODB, it reconfigures the trigger electronics and operates the 
beam blocker (for cosmic runs) automatically. At the end of a calibration run his- 
tograms are automatically evaluated to obtain pedestal values, software gains for 
ADC and timing offsets for TDC data. 



12.2 Single event display 



The PIBETA single event display (SED) was written in the Java programming 
language. It has proven to be a valuable debugging tool for verifying the correct 
operation of the user analyzer code. It has also proven to be indispensable in the 
development of the DSC signal fitting algorithms. The display program can run 
online by reading from the data stream as well as offline by taking the input from 
raw MIDAS files. The graphics interface menu offers the following display choices: 

(1) a segmented Csl calorimeter map in the Mercator projection; 

(2) the front or side cross-sections of the detector with the active target, MWPCs, 
plastic veto and the Csl calorimeter; 

(3) a digitized waveform from a single DSC channel. 

The calibrated ADC/TDC values, i.e., energies and event times, can be shown asso- 
ciated with the individual detector shapes, like PV segments and Csl modules. The 
cross-sectional views show the reconstructed charged particle tracks, while different 
particle types (positrons, protons, photons and cosmics) are identified by differently 
colored line types. 
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The DSC option shows both the digitized waveform of the selected channel and the 
sum of fitted system functions in a 125 ns time window. The number of hits found 
is indicated by circles located at the DSC TDC timings. 

Program options include the selection of a particular trigger or an arbitrary trigger 
combination. The scanned events can be displayed sequentially, or, alternatively, 
triggers can be stepped through in a loop. A single selected event can also be 
specified by its serial number. 

The top panel of Fig. 39 shows the SED of a /i"*" e^U(>Ufj_e^ e~ event in a lateral 
{x-y) cross section through the detector. The bottom panel illustrates a digitized 
waveform of one active target segment (ch. 283) displayed in a SED window. In 
this example the fitting routine finds two hits, one associated with a stopping pion 
signal (first, larger peak in channel 24, i.e., at -48 ns), followed by a stopping muon 
pulse from the tt'^ — > n'^Vfj, decay (second pulse, in channel 42, i.e., at •84ns). 



12.3 Experiment WWW control page 



Detector operation is monitored remotely at all times by experiment collaborators 
who could be residing in different world time zones. Monitoring is facilitated by 
the implementation of the WWW interface to the online experiment. At least two 
shift takers were usually present locally at PSI during production running. Off-site 
workers followed the experiment via DSL or modem lines. 

At a glance, the WWW experiment control page provides an overview of status of 
the experiment, as well as the conditions of the detector subsystems. The data taking 
can be started, stopped, paused or resumed from the control page. The experiment 
URL is linked to the ODB and CAMAC CNAF command page, a window with the 
latest automatically generated messages from the experiment tasks, the electronic 
logbook system, alarm system page, the menu of all experiment tasks, the history 
page and the online help link. 

A separate menu line leads an experimentalist to the trigger settings, detector rates 
and ratios of rates, the PSI accelerator status display, the WebPAW interface, and 
the online PIBETA detector handbook. A user can also examine trigger and scaler 
variables, check all demand and measured high voltages, as well as the wire chambers 
and the area beam line settings. Also highlighted, due to their importance, are 
the activities of the logger channels, the sizes of data files for the current run, 
the available free disk space and the progress in writing the data to the magnetic 
tape and the remote ftp archive. About half of all detector malfunctions can be 
successfully corrected by user intervention from a remote location. 
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12-4 Electronic logbook 



A traditional paper logbook was not a practical option for a number of reasons. To 
adress this, the MIDAS system contains an electronic logbook. Shift workers can 
enter logbook messages with optional attachments through the same Web interface 
used for experiment control. Online histograms can be attached in addition to the 
variables defined in the online database. Messages are classified into certain groups, 
which makes it simple to search for a specific entry. In addition, the analyzer has 
access to the electronic logbook and can submit regular entries concerning, for 
example, the automated gain calibration. 

Use of the electronic logbook greatly enhanced access to the experiment and mutual 
communication among collaborators both on-site and remotely. The e-logbook files 
were copied on a daily basis between PSI and the University of Virginia for local 
access and backup. A stand-alone version of the electronic logbook is available for 
offline browsing of the logbook entries [45]. 



12.5 Alarm system 



To optimize the efficiency of data taking, malfunctions of the detector or the beam 
line must be quickly detected and promptly remedied. MIDAS contains an alarm 
system for that purpose. This system regularly checks the operation of all experi- 
ment computers and programs. In addition, any parameter value defined by the slow 
control system and/or known to the analyzer can be verified. We use this facility to 
check: 

(1) temperatures in the detector enclosure and the electronics hut; 

(2) gas pressure in the MWPC system; 

(3) detection efficiency of the MWPCs; 

(4) operation of the tape drives and uninterrupted archive backup; 

(5) rates of the pion beam stops and triggers; 

(6) scaler rates for all Csl crystals and plastic veto counters. 

These computerized checks ensure good stability of the experiment over long mea- 
surement periods. The detector rates are especially sensitive to beam line problems 
and electronic unit failures. In case of an alarm, the shift persons, both local and 
remote, are notified by telephone calls through a modem connected to the back- 
end computer. This arrangement makes it possible to run the detector unattended 
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over certain periods while retaining the abihty to react to problems in less than 15 
minutes at any time. 



13 Energy calibration of the PIBETA calorimeter 



Energy calibration of the PIBETA calorimeter involved two correlated processes: 
equalizing discriminator thresholds (expressed in MeV) for 220 Csl detector signals 
that define the calorimeter trigger and calibrating signal gains of 240 Csl detec- 
tors at the opposite analog ADC/DSC branch. We accomplished first the trigger 
branch threshold equalization by adjusting the high voltages of the individual Csl 
PMTs. Once all Csl discriminator thresholds, expressed in calibrated MeV units, 
were matched (usually within ±1 V, or 0.5 % in signal gain), the Csl analog signals 
measured in the FASTBUS ADC branch were matched by introducing software 
gain factors for individual detectors. These two adjustments can be done both man- 
ually by the experiment operators on duty, as well as automatically by a computer 
program. 



13.1 Calorimeter PMT gain adjustment 



Using a digital oscilloscope we confirmed that the individual daisy-chained UVA 125 
cluster discriminator /summers all have identical demand threshold levels (within 
1 %). The difference (if any) between a signal gain of an individual Csl detector at 
the trigger branch on one side, and at the input to a digitizing ADC/DSC system on 
the other side, could arise from different signal attenuations in the delay cables of two 
branches, differences in the resistor values in UVA 126 and UVA 139 passive signal 
splitters and/or small mis-matches in time offsets of the two branches. We measured 
ratios of the fast trigger pulses to delayed pulse charges using an oscilloscope in the 
single shot mode and found the normalized ratios to vary up to 10 % for the same 
Csl detector. The conclusion is that once Csl signals are perfectly gain-matched at 
the trigger inputs, thus giving the best overall trigger energy and timing resolution, 
the signals in the ADC branch have to be simultaneously scaled by 0.9-1.1 to attain 
the best possible calorimeter energy resolution. 

For a detector equipped with an n-stage PMT, a normalized real gain change g, 
relating two different configurations 1 and 2, depends on the ratio of the software 
gains Si and the corresponding high voltages HVj: 




(20) 
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The above equation relates the net gain factor of a scintillator detector in configu- 
ration 2 to the gain in configuration 1. 

To facilitate and automate the detector energy calibration, an analyzer routine was 
developed that could optionally: 

(1) apply an overall gain factor by varying all Csl PMT high voltages or software 
gain multipliers; 

(2) adjust the Csl ADC gains at the trigger discriminators via HV changes by 
comparing the counting rates in the tt — > ez^ peak and the high energy edge of 
the /J, — > evv spectrum; 

(3) match the Csl ADC gains using the fitted tt — > ei/ peak histograms; 

(4) match the Csl ADC gains using the fitted n — > evv energy spectrum; 

(4) adjust the PV detector geometric-mean ADC gains using the fitted peaks of 
minimum ionizing e^s. 

At detector start-up time, the gain-match routine is run manually by an experienced 
operator. Several iterations of the Csl PMT high voltage adjustments are necessary 
before the /j,'^ — > e'^UeUn to 7r+ e+i/g counting rates for the high threshold trigger 
are fixed to the default 3:1 ratio. The software gains of individual Csl detectors are 
then set by matching the positions of the 69.8 MeV tt"*" —>■ e^Ve peaks. Consequently, 
all Csl high voltages are adjusted automatically by the computer program when the 
preset cumulative statistics in the reference histograms are attained. This dynamic 
gain monitoring occurs daily and the gain-matching routine automatically submits 
the appropriate e-log entries. 



13.2 Calorimeter energy resolution 



We calibrated the 5-150 MeV dynamic range of the calorimeter by employing 

(1) monoenergetic 69.8 MeV 7r+ e+t'e positrons; 

(2) the energy line shape of reconstructed 7r°s following the tt/? decay at rest; 

(3) the energy line shape of the SCX tt^s produced on the CH2 target; 

(4) the jU"^ — ^ e^Ve^n energy spectrum; 

(5) the cosmic muon energy spectrum. 
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The most important calibration point comes from the vr"*" e^v^ line: the exper- 
imental spectrum shown in Fig. 40 is used to fix the single conversion factor to 
the absolute energy scale as well as to fine-tune the parameters (mostly the photo- 
electron statistics) in the detector Monte Carlo simulation. Once established, these 
parameters are "frozen" in the offline data analysis and MC simulation. The energy 
distributions for all other reaction channels and probes are then required to agree 
independently with the same GEANT MC simulation. 

The tt"*" — > e^Uf, calorimeter peak position is determined by energy losses in the 
active target, plastic veto scintillator, and the insensitive layers in front of the 
Csl crystals, positron annihilation losses, photoelectron statistics of individual Csl 
modules, and axial and transverse coefficients parameterizing the nonuniformities of 
Csl light collection. At the measured peak position of 62.55 ± 0.03 MeV, the FWHM 
fractional resolution AE/E is 12.8 ± 0.1 %. 

The demonstration that the energy response of the calorimeter is simulated correctly 
also at the lower energy scale between 5 MeV and 50 MeV is given by the Michel 
/i"'" e^Uf>i>fj_ energy spectrum. The experimental histogram accumulated for 10 
runs as well as the expected Monte Carlo spectrum are shown in Fig. 41. The 
per degree of freedom is 1.3. 

The precise trigger-defining absolute high energy and low energy discriminator 
thresholds can also be determined from the Michel energy spectra. We find that 
the average values in the production running were £'ht = 51.8 it 0.1 MeV and 
EisT = 4.5 lb 0.1. The corresponding threshold widths, important in the correct sim- 
ulation of the detector response were measured to be ctht = 2.4 ± 0.1 MeV, and 
(Tlx = 1-2 ± 0.1 MeV, respectively. 

Fig. 42 shows the experimental vr" spectra. The superimposed GEANT simulation of 
the tt" line shapes agrees very well with the data. The calorimeter fractional FWHM 
resolution at the energy of 128.70 ± 0.06 MeV is l^E/E = 10.0 ± 0.1 %. 



14 Timing response of PIBETA calorimeter 



The calorimeter time resolution depends on the intrinsic time resolution of the 
individual Csl modules, the spread in the arrival of analog PMT signals at the 
trigger point where the analog Csl summing is done, and the uncertainties of the 
software time offsets. Before assembling the calorimeter we measured the intrinsic 
time resolutions of all component Csl modules using cosmic muons as a probe. Csl 
times are determined relative to a small plastic scintillator counter. The average 
Csl detector rms TDC resolution specified in such a way is 0.68 ns. The details of 
these measurements are provided in Ref. [31]. 
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14.1 Multi-hit TDC data 



The discriminated signals of all PMTs are timed in 96-channel FASTBUS LeCroy 
Model 1877 multi-hit TDC units. The instruments were operated in common stop 
mode with 0.5 ns least significant count, measuring both leading and trailing edge 
time information over an interval of 512 ns. Up to 10 sequential leading and trailing 
edge hits per event per TDC channel are stored in the analyzer bank and, optionally, 
in the corresponding column-wise N-tuple. However, the physics analysis dictates 
a more convenient way to present these data. We chose a scheme which saves in- 
formation on up to three hits in any channel (leading edge time plus width) along 
with the information on the number of hits in the channel. The scheme is flexible 
and can easily be modified to record more information. 



14-2 Trigger timing 



The timing spread of the 220 trigger-defining Csl analog signals is checked period- 
ically in timing calibration runs with the prompt trigger. The idea is to find the 
time difference between a single reference detector, in our case the active degrader, 
and each Csl counter. This type of timing histogram, associated with a given Csl 
detector, is incremented only if a charged particle track is identified as a fast pro- 
ton (Ep > 60 MeV) in the plastic veto hodoscope and 80 % of the shower energy is 
contained in that module. The total energy contained by the calorimeter is used to 
calculate the time-of-flight correction, a term that is as large as 1.0 ns for 100 MeV 
protons. We use the proton events because ^ 1 % statistics in TDC spectra is ac- 
quired within one hour of data taking. The peaks of the timing histograms are fitted 
at the end of the run and the peak positions are ordered relative to the slowest Csl 
detector. The resulting information is used to add trigger cable delays, available in 
0.5 ns increments, to the faster Csl lines, at the patch panel preceding Csl discrimi- 
nator/summer units. Three iterations of this procedure resulted in a 0.86 ns relative 
trigger rms timing spread (Fig. 43). 



14-3 TDC calibration: zero offsets and slewing 



TDC calibration is accomplished via two independent corrections, both applied in 
software. The primary TDC correction compensates for different cable delays of the 
digitizing branch. The zero time is defined as the center of gravity of the self-timing 
peak for each detector channel. The self-timing peaks are evaluated in the end-of- 
run routine of a special analyzer module. If the appropriate analyzer flag in the 
ODB is enabled, all raw TDC histograms are fitted with Gaussian functions. The 
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results, namely timing peak offsets and the rms values of the peaks, are stored in 
arrays in the ODB database. Since the algorithm needs a clean self-timing peak, this 
procedure is performed only on the prompt event histograms. All raw TDC values 
are then modified in the analyzer by subtracting the zero point offset. Negative 
TDC values clearly correspond to signals which come before the trigger, positive 
time values are associated with the delayed calorimeter events. 

The secondary TDC correction linearizes the slewing of TDC time caused by the 
differing amplitudes of ADC signals. A smaller amplitude signal takes more time to 
rise to the fixed discriminator threshold than a larger signal. The result is an artificial 
energy dependence of TDC values with lower energy signals registering later times. 
The data for each individual channel are fitted to determine the TDC versus ADC 
dependence. The secondary TDC correction is implemented in offline analysis by 
subtracting the fitted function from each TDC value. We have implemented this 
solution using the basic functional form: 

CTDC = TDCo + a ■ (ADC - b)" (21) 

where TDCq, a, b, and c are free parameters of the fit, ADC is the calibrated ADC 
value proportional to the deposited energy, and the result, CTDC is the corrected 
TDC value. Fig. 44 shows the energy dependence of one representative Csl TDC 
and the reduction in the time slewing after applying the correction. This procedure 
is automated in the analyzer and includes all PMT signals. The experimental data 
are fitted automatically and the parameters from the above equation are stored in 
the ODB database. 

We use the coincident photon pairs following the decay of tt^s copiously produced 
by single charge exchange (SCX) inside the degrader and target material to quan- 
tify the improvement in the time resolution due to the algorithm described above. 
Fig. 45 shows the relative timing between two vr^ photons before and after the 
correction. Both distributions are fitted with the sum of a constant term and a 
Gaussian function. The algorithm reduces the rms width from 1.30 ns to 0.84 ns, 
or 0.59 ns per single photon shower. The signal-to-background ratio, defined as the 
ratio of the height of the Gaussian peak P2 to the constant term PI, increases more 
than fiftyfold, from 24 to 1214. 

Fig. 46 shows the dependence of the online timing resolution on the 7r+ stopping rate. 
The top panel is the radiative vr"*" decay signal recorded with one-arm calorimeter 
trigger at production rate of ~ 7.3 • IO^tt^/s. The same e+-7 timing difference 
shown on the bottom panel corresponds to the data for the lower intensity beam 
of ~ 1.4 • 10^7r+/s. The low-beam rms calorimeter timing resolution is again 0.6 ns 
per shower. 

We have continuously monitored the absolute t = point using the precise accelera- 
tor primary beam structure. The precise PSI accelerator frequency is 50.63280(4) MHz, 
and the phase stability of the primary quartz oscillator is better than 0.01°. This 
stability translates to the time interval between pulses of 19.750 ns with ATj-f = 
0.0028 %. The proton beam bunch rms width is ~ 0.03 ns. 
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We have studied the timing offset accuracy that we can achieve relying on SCX tt^ 
event timings. The relevant variable defining the timing of SCX vr^'s relative to a 
Tr~^ stop signal is: 

[m,+^,)/2]-t^G, (22) 



where and t^^ arc the TDC values for two neutral Csl calorimeter showers that 
reconstruct to vr" and vrG is a TDC datum of the vr"*" gate signal. The gate width 
was set to ~ 190 ns and thus captures 10 beam pulses. The beam pulse closest to 
t = is suppressed in the trigger by the beam veto; therefore we do not use the 
pion stop peak in the analysis. 

As an example of our approach we show the SCX 7r° timing peaks for a set of 236 
runs in Fig. 47. We fit the TDC lineshapes with Gaussian functions and extract 
the peak positions and widths. The average non-linearized timing peak width is 
a = 1.493 ± 0.014 ns. A least-squares MINUIT fit of the extracted peak times of the 
9 beam bursts with a period Trf left as a free parameter gives: 

to = -4.003 ± 0.026, and = 19.751 ± 0.010. (23) 



15 Angulcir resolution of PIBETA calorimeter 



The calorimeter granularity is such that if a module is hit centrally, it will contain 
about 90 % of the shower energy. Three Csl modules at most should receive signifi- 
cant portions of a single particle's shower energy. The median number of calorimeter 
modules with energy depositions above 0.25 MeV recorded in low beam flux runs 
(with no accidental pile-up) with a one-arm high threshold trigger is 25. 

The angular direction {9c, (j)c) of a particle initiating an electromagnetic shower in 
the PIBETA calorimeter is found as an energy-weighted mean: 

N 

Oc = ^ , (24) 

i=0 



where 6i is the polar angle of an individual Csl module, and N is the number of 
nearest neighbors of the crystal with maximum energy deposition Ef^"-^. The sum 
is over the central crystal "0" and all of its nearest neighbors and thus contains 6 
to 8 terms, depending on the shape of the centrally hit module. The formula for 
the azimuthal angle (pc is identical to that for the polar angle with the replacement 
Oi (t)i. 
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Three different weighting functions LOmi{Ei) are examined (see Refs. [46-54]): 



where a and ao are empirical parameters and Etot is the total energy deposited in 
the N + 1 modules that define the shower cluster. The linear weighting formula 
predicts an rms angular resolution of 2.2° when coded into a GEANT Monte Carlo 
that generates 70MeV positrons emanating from the target center. A slightly bet- 
ter resolution of 2.0° is achieved by power weighting with a = 0.7, though both 
algorithms overemphasize the calorimeter granularity. The logarithmic weighting 
method of shower localization [52] is motivated by the exponential fall-off of the 
shower's transverse energy deposition [50]. The log- weighting algorithm gives the 
best directional resolution and it was employed in the offline analyzer code. 

The optimum values for the power-weighting exponent a and the logarithmic al- 
gorithm parameter ao are found by two different methods: GEANT Monte Carlo 
simulations of the detector uniformly illuminated by positrons and photons in the 
energy range 30-80 MeV and comparison of the particle directions reconstructed 
from the calorimeter response to 70 MeV 7re2 positrons and the positron MWPC 
tracks. The GEANT optimization gives a = 0.5 and oq = 5.5 for the best overall an- 
gular resolution. Fig. 48 shows the distribution of differences between the charged 
particle direction deduced from MWPC hits and the direction found from power- 
weighted calorimeter energies. The rms value of 2.0° between the MWPC tracks and 
calorimeter shower centroids agrees well with the Monte Carlo prediction quoted 
above. 



16 Energy clustering and tracking algorithm 



The cluster finding algorithm operates in the online analyzer program and makes 
the first attempt at identifying calorimeter shower clumps due to the interaction of a 
single particle. The initial inspiration for our algorithm was the Crystal Ball idea of 
"a clump discriminator function" [28]. The algorithm first constructs a list of up to 
6 nonadjacent calorimeter modules in order of decreasing deposited energy ("clump 
centers"). The minimum calibrated energy allowed for a clump center in online 
running is 1.0 MeV, which is lower than the low energy discriminator threshold of 
~4.5 MeV. 

The clump energies Ed are calculated next by adding the calibrated energies of 
neighboring modules to the energy recorded in the clump center, provided that the 
TDC hits of neighbors fall within a specified time window. The default value of 
the window width is ±5 ns with respect to a clump center time. Energies of crystals 
outside that time window are not included in the energy sum, under the assumption 




LOli = Ei linear weighting, 

uj2i = Ef power weighting, 

Wii = max[0, ao + ln(£^i) — In^Etot)] logarithmic weighting. 



(25) 
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that they arc related to accidental coincidences. The number of nearest neighbors 
varies from 5 for a centrally hit pentagon, to 7 neighbors for the outlying HEX-C 
modules. When acquiring data at lower beam intensities, with correspondingly lower 
accidental pile-up rates, a provision is made for two rings of neighboring modules 
to be included in the clump sum. The time associated with a clump is calculated as 
the energy weighted average of all clump members. The clumping algorithm finally 
saves energies, times and the relative angles between all clump pairs in an analyzer 
bank. A similar bank is also created for up to 6 plastic veto "clumps" which contain 
energies and times of single hit PV detectors. 

The tracking algorithm begins by finding the "best" pairs of hits in the inner and 
outer MWPC chambers, using the criterion that the most probable tracks point to 
the stopping target in the center of the detector. The charged tracks identified in this 
way are then associated with the Csl and PV clumps defined previously. Finally, all 
Csl clumps that do not pair with the charged MWPC tracks are defined as neutral 
tracks. The single track information sets the stage for particle identification. The 
target, PV and Csl clump energies (Eqs. 2 and 3) and corresponding times, and 
MWPC directions of the identified charged tracks (if any) are used to associate a 
positron, proton, photon, neutron or cosmic muon label with each track. This final 
piece of information is appended to the track analyzer bank. 



17 Detector Monte Carlo description 



A complete GEANT Monte Carlo description of the PIBETA detector that includes 

all major sensitive as well as passive detector components was developed [11]. The 
user code is written in modular form in standard FORTRAN 77 and organized into 
over 300 subroutines and data files [55] . A graphical user interface to the code based 
on Tel routines supplements the program. 

Standard GEANT 3.21 routines are used to describe the electromagnetic shower pro- 
cesses in Csl (leak-through, lateral spreading, and backsplash) as well as the electro- 
magnetic interactions in the active target and tracking detectors (bremsstrahlung, 
Bhabha scattering, and in-flight annihilation for the positrons, and pair production, 
Compton scattering, etc. for photons). 

The individual PIBETA detector components, both passive and active, can be repo- 
sitioned or disabled without recompiling the code. The following detectors are de- 
fined in the Monte Carlo geometry: 

(i) a beam counter, active plastic degrader and several versions of segmented 
stopping targets; 

(ii) two concentric cylindrical multiwire proportional chambers used for charged 
particle tracking; 
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(iii) a 20-piece cylindrical plastic hodoscope used for charged particle discrimina- 
tion; 

(iv) a 240-module pure Csl calorimeter sphere; 

(v) a 5-plate cosmic muon veto scintillator system; 

(vi) a passive calorimeter stand and individual detector support systems; 

(vii) detector phototubes and HV divider bases, and a lead brick shielding. 

User input to the simulation code requires the detector version, e.g., the run year, 
beam properties, selected reaction/decay final states, and the version of ADC and 
TDC simulation codes. A user can select one particular final state or any combina- 
tion of different final states with the relative probabilities defined by the reaction 
cross sections and decay branching ratios. The simulation of ADC values and TDC 
hits accounts for individual detector photoelectron statistics, axial and transverse 
light collection nonuniformities, ADC pedestal variations, electronics noise, and 
event pile-up eff'ects. Photoelectron statistics and individual detector light collection 
nonuniforniity coefficients are initialized from the RASTA database file ([31]). The 
calculation of accidental coincidences up to fourth order is optional. 

Selectable options include the effects of: 

(1) lateral and axial extent and the divergences of the stopping pion beam; 

(2) positron and muon beam contamination; 

(3) photonuclear and electron knockout reactions in Csl material [56]; 

(4) aluminized Mylar wrapping of Csl modules and plastic veto staves; 

(5) the gaps between the Csl detector modules; 

(6) the temperature-dependent light output coefficients of the individual calorime- 
ter modules; 

(7) gain instability and drifts of the calorimeter detector modules; 

(8) electronic discriminator thresholds and ADC gate widths. 

Particular attention was paid to the correct accounting of individual Csl module 

software gains. The set of gain constants depends on the optical properties of in- 
dividual Csl crystals as defined in our GEANT database, but also upon the incident 
particle chosen for the Monte Carlo calibration runs and its energy because of differ- 
ences in shower developments of photons and positrons. In the simulation calculation 
the values of the detector software gains allow the same user control as the detector 
high voltages in the operation of the physical detector. The software gains were 
determined in an iterative procedure constrained by the Monte Carlo tt"*" — > e^u^ 
positron ADC spectrum in each Csl detector. Similarly, in the real experiment, the 
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high voltages of the individual Csl detectors are also set by matching the positions 
of the 69.8 MeV positron peaks. 

The output of the calculation saves the selected experimental layout and the fi- 
nal states chosen in the calculation, as well as the production cross sections and 
branching ratios used. A number of physical variables that could be of interest are 
optionally filled in 1- and 2-dimensional histograms. Row-wise and Column-wise PAW 
N-tuples [42] are used to digitize individual events. Simulated energy depositions 
and the ADC and TDC values associated with the sensitive detectors are saved on 
an cvcnt-by-cvcnt basis. The single event display operating on the output can be 
used to examine a particular event in detail. 



18 Radiation hardness and aging of active detector elements 



During production running in 1999 and 2000 the detector operated in the pion beam 
for a total of 297 days accumulating radiation doses of up to 2 • 10^ rads. During 
that period the active target registered a total of 1.4 ■ 10^^ tt"*" stops. Table 5 lists 
the radiation doses received by individual active detector components. The more 
detailed analysis of the radiation resistance of the PIBETA detector is presented in 
a separate publication [57]. 

The gain of the forward beam counter decreased approximately linearly over a period 
of one calendar year due to exposure to 2 Mrad of radiation. The decrease in gain 
was ~20%/Mrad. The energy resolution for through-going pions was also slightly 
degraded, changing from 13.1% to 13.9%. 

In the course of data taking the active degrader counter was irradiated with a 
1.4 Mrad dose, producing an average gain decrease of 15%/Mrad. The temporal 
degradation in the counter energy resolution is much smaller but is still noticeable. 

The five central target segments AT0-AT4 received the bulk of the beam exposure, 
and consequently show the largest temporal variation in the gain factor and energy 
resolution over the 15 month period. Table 5 details the responses of three different 
segments. The cumulative absorbed radiation dose is calculated to be 0.5 Mrad per 
segment. The reduction in the gain of the central target scintillator is considerable 
(~52%/Mrad). The energy resolution was degraded by a full 80%. This level of 
radiation damage forced us to manufacture two additional copies of the active target 
detector. The target assembly was replaced on an annual basis. 

The temporal dependence of the light output was similar for all 20 PV hodoscope 
staves. Table 5 lists the data for two individual PV detectors as well as the average 
numbers for the complete PV system. The gain factor dropped by 0.25%/krad, 
giving a cumulative degradation over the production period of 10%. The energy 
resolution on the other hand shows a barely measurable change. 
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The overall Csl calorimeter gain and energy resolution were determined from cal- 
ibrated energy spectra of 240 Csl detectors. The results are again summarized in 
Tables 5. The PIBETA calorimeter gain decreased 17% and the overall online rms 
energy resolution for the 70MeV e+s was lowered from 5.5% to 6.0%. For an illus- 
tration we present the initial and final gain and resolution values for 6 individual 
Csl detectors. 

In conclusion, we find measurable decreases in energy gains and degradations in 
energy resolutions of the PIBETA detector subsystems, in particular for the beam 
detectors that were exposed to high radiation doses as well as for the more radiation- 
sensitive pure Csl calorimeter modules. These changes in detector responses are 
monitored online throughout the production running and documented in the replay 
data analysis. The changes affect the energy calibration of the PIBETA detector 
elements and have to be taken into account properly when defining the energy cuts 
in the physics analysis of rare pion and muon decays. 



19 Results and conclusions 



We designed, built and commissioned the PIBETA detector — a nonmagnetic large 
solid angle detector with an active stopping target. It was used for studies of rare 
pion and muon decays at PSI. In this paper we have provided a comprehensive 
description of the design of the individual detector subsystems and their calibration 
and performance in the intense tt^ beam. 

The major performance parameters of the detector's Csl calorimeter are: 

(1) the rms energy resolution for TOMeV positrons i^EjE of ~5.0%; 

(2) the rms energy resolution for decay-at-rest 7r°s l^EjE of ~4.0%; 

(3) the two-arm relative time resolution for two 7r°-decay-at-rest photons ~ 1.3 ns 
online and ~ 0.8 ns offline; per-shower timing resolution ~ 0.6 ns; 

(4) the angular resolution of TOMeV positron tracks ~2.0°. 

The detector's charged particle tracking system incorporating a pair of cylindrical 
MWPCs and a plastic veto hodoscope has a stable charged particle detection inef- 
ficiency of < 2 • 10~^. Positron identification, when used in conjunction with data 
from Csl calorimeter, is better than 99.8% accurate. 

The pion beam geometry and beam transmission, temperature and humidity inside 
the detector housing, relative trigger rates, detector time offsets and PMT gain 
factors are proven to be very stable during the three calendar years of the data 
acquisition. Detector operation is nearly 100% automated, requiring only a few 
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experimenters to be physically present and on call at the PSI site during the routine 
data acquisition, with other collaborators working remotely. 

The apparatus can easily handle pion stopping rates of IMHz and lower. It can 
detect virtually background free rare tt(3 decay events, as well as radiative pion 
and muon decays. Minimal backgrounds at the highest stopping rate arise from 
accidental coincidences with ordinary muon decay events and cosmic muons. So far 
we have collected physics and calibration data based on 2.3 • 10^^ stopping pions. 
The analysis of the vr/? data sample and radiative pion and muon decay events is 
underway and will soon be ready for publication. 
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Fig. 1. A schematic cross section of the PIBETA apparatus showing the main com- 
ponents: forward beam counter (BC), two active colHmators (ACl, AC2), active 
degrader (AD), active target (AT), two MWPCs and its support, plastic veto de- 
tectors (PV) and PMTs, pure Csl calorimeter and PMTs. 

Fig. 2. A photograph of the PIBETA assembly located in the irEl experimental 
area: the air-conditioned enclosure housing the PIBETA calorimeter is in the center, 
focusing quadrupole magnets are on the left, and the lead-wall cosmic house is 
opened at right. The analog delay cables, high voltage supplies and air-conditioning 
system are partially blocked from view by the cement shielding at front, while the 
fast electronics house is visible at the back of the photo. 

Fig. 3. A layout of the ttEI secondary beam line in the main experimental hall 
of the Paul Scherrer Institute's ring accelerator. Focusing quadrupole magnets are 
labeled with "Q" , bending dipoles are labeled with acronyms beginning with "A" , 
and passive collimators are labeled with acronyms beginning with "F" . 

Fig. 4. A layout of the PSI ttEI experimental area for the PIBETA experiment (top 
view). The PIBETA detector is located on the platform which also carries the fast 
electronics, air conditioning system, high voltage supplies, and analog cable delays. 
The focusing quadrupole magnets in the beam line channel (Qs) , the forward beam 
counter (BC) and the passive collimator (PC) are also indicated. 

Fig. 5. The momentum spectrum of the positive pions at the face of the active 
degrader as calculated by the TURTLE beam transport program. 

Fig. 6. The pion signal ("7r-stop") is defined as a four-fold coincidence between a 
forward beam counter (channel 1), an active degrader (2), an active target (3), 
and a 19.75 ns rf cyclotron signal (4). In the above digital oscilloscope snapshot the 
channel 4 scale shows three rf periods. 

Fig. 7. A Polaroid photograph of the beam spot taken in front of the forward passive 
collimator PC. The brighter ellipsoid right of the center is the e"*" beam stopped in 
the collimator material. The fainter spot in the center is the momentum-analyzed 
tt"*" beam. The pions and positrons at that point are separated by 40 mm in the 
horizontal plane. 

Fig. 8. A schematic drawing of the forward beam counter inside the light-tight 
housing. A feedthrough and a single photomultiplier tube are drawn only on the 
right side. The beam direction is perpendicular to the plane of paper. 

Fig. 9. The uncalibrated pulse-height spectrum of the 40.6 McV incident vr^ beam 
as detected in the thin forward beam counter. Gaussian fit parameters for the lower 
energy part of the 7r+ peak are showed in the statistics window. 
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Fig. 10. A cross sectional view (left) and a frontal view (right) of the active degrader 
counter. The degrader counter has the shape of a beaker vessel. Beam pions incident 
parallel to the central detector axis arc moderated in 30 mm of scintillator material 
before stopping in the center of active target. The four light guides which end up 
in a cross section of 6 x 6 mm^ are also shown. 

Fig. 11. The uncalibrated ADC spectrum of moderated vr+s in the active degrader 
scintillator. A 25 ns ADC gate is used with a 10-bit CAMAC unit. The Gaussian fit 
is superimposed on the data points. 

Fig. 12. A sketch of the regular PIBETA active target, composed of 9 detectors. 
The segment sizes are chosen to balance the scaler counting rate. 

Fig. 13. A photograph of the regular PIBETA active target. Acrylic tapered light 
guides are glued to 9 optically isolated target segments. 

Fig. 14. The pulse-height spectrum of pions stopping in the central active target 
segment. The peak-to-tail area ratio depends on the beam divergence. 

Fig. 15. Reconstructed 2-dimensional shape of the 7r+ beam superimposed on an 
outline of the segmented active target. Data for target counting rates are collected 
during one 48 hour long series of runs. 

Fig. 16. The longitudinal (z) component of the tt"^ stopping distribution in the active 
target calculated in a GEANT simulation using the momentum distribution given in 
Fig. 5. 

Fig. 17. Beam particle identification using the time-of-flight (TOF) method. The 
TOF difference between the forward beam counter BC and the active target AT, 
which are separated by 3.87 m is plotted along the vertical axis. The horizontal axis 
shows the particle TOF between the production target E tagged by the 50 MHz 
accelerator rf signal and the active target AT. The beam path length between these 
two points is 16.83 m. 

Fig. 18. MWPCi angular resolutions in the azimuthal angle (pcH (top panel), and 
axial coordinate zqh (bottom panel) extracted from a sample of cosmic muon events. 
The data points for Acj) and Az resolutions are fitted with a sum of two Gaussians, 
the broader one accounting for the cosmic muon scattering in the apparatus.. 

Fig. 19. The MWPCi detection efficiency for minimum ionizing charged particles 
as a function of the azimuthal angle (top panel). The bottom panel shows the 
equivalent histogram for the outer chamber MWPC2. 
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Fig. 20. A GEANT rendition of the plastic veto hodoscope array surrounding the 

active target in the center of the PIBETA detector (top panel). The dimensioned 
side and end views of one plastic veto scintillator stave with the light guides glued 
at the both ends are shown in the middle and bottom panels, respectively. 

Fig. 21. The calibrated energy spectrum of minimum ionizing positrons (left peak) 
and 5-150 MeV protons (right peak) in the PIBETA plastic veto hodoscope. 

Fig. 22. The AE-E particle discrimination between protons and minimum ionizing 
positrons using the calibrated ADC values of plastic veto and Csl calorimeter. 

Fig. 23. The detection efficiency of the PIBETA plastic veto hodoscope for the 
minimum ionizing charged particles plotted as a function of the azimuthal angle. 

Fig. 24. An expanded view of the detection efficiency of the plastic veto hodoscope as 
a function of the azimuthal angle: the inefficient regions correspond to the ~0.08 mm 
gaps between the 20 individual PV staves. 

Fig. 25. Nine different Csl scintillator shapes that make up the PIBETA calorimeter: 
(i) PENTA, (u) HEX-A, (in) HEX-B, (iv) HEX-C, (v) HEX-D, (vi) HEXDl, (vu) 
HEXD2, (viii) VETOl, and (ix) VET02. The calorimeter proper is built from 220 
Csl crystals while 20 Csl veto crystals surround the beam entry and exit regions. 

Fig. 26. The top panel shows the geometry of the pure Csl shower calorimeter. The 

sphere is made up of 240 elements, truncated hexagonal, pentagonal, and trapezoidal 
pyramids; it covers about 77% of Ivr in solid angle. The bottom panel separates 
one supcrcluster and its complement. The PIBETA calorimeter is made out of 10 
nonover lapping super clusters. 

Fig. 27. The breakdown of the PIBETA calorimeter in the Mercator projection with 
the numbering scheme adopted for the detector clusters and superclusters. The top 
indices refer to the module shapes, the bottom labels represent the crystal serial 
numbers. The beam orientation as well as the polar and azimuthal coordinate axes 
are indicated. 

Fig. 28. Two drawings of the cosmic veto counters: the top panel shows the make-up 
of a single scintillator plane. The bottom panel shows the details of the light readout 
geometry via wavelength-shifting bars. 

Fig. 29. The energy and time resolution of the cosmic muon veto detector. The CV 
timing with respect to a calorimeter TDC hit is plotted against the uncalibrated 
CV ADC values. 

Fig. 30. (i) Cosmic muon low-threshold trigger: the energy spectrum of the most 
energetic clump [Ei > 5 MeV) in the Csl calorimeter, (ii) Energy sum of the two 
high threshold Csl clumps (£^ci,C2 > 52 MeV) with a relative opening angle greater 
than 150° (log scale). The Monte Carlo predictions are shown with full line functions. 
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Fig. 31. The temperature stability of the pure Csl calorimeter: temperature sensor 
readings in 10 seconds intervals spanning 100 hours are shown. The maximum tem- 
perature variation of <0.1° C corresponds to a variation of <0.15% in the Csl light 
output. 

Fig. 32. "The system function" describing a default vr"*" waveform in the active 
degrader counter. The functions of that type are derived by averaging 100,000 single 
hit waveforms in the 250 ns DSC gate. 

Fig. 33. Comparison between the DSC-derived data and the conventional ADC and 
TDC measurements. Detailed interpretation is given in the text. 

Fig. 34. The calibrated ADC spectrum of a single Csl detector measured using a 

random event trigger at the beam intensity of 8 • lO^vr^/s. The mean accidental 
energy deposited in a single Csl detector averaged over all calorimeter modules at 
that beam flux is 0.15 MeV. 

Fig. 35. A block diagram of the fast trigger electronics. The PIBETA trigger logic 
is implemented in three PSI LB500 logic units that use the "on-board" repro- 
grammable logic device ispLSI 2128E-LT165 from Lattice Semiconductor [40]. 

Fig. 36. A digital oscilloscope snap-shot of a vr-stop signal (channel 1), a pion gate 
(2), a beam signal (3), and a "Csl high threshold" signal (4). The "tt — >■ ev high" 
calorimeter trigger is defined as a 3-fold coincidence 2 • 3 • 4. 

Fig. 37. The calibrated pulse-height spectrum of the most energetic Csl "clump" 
for the events identified in the plastic veto hodoscope as prompt protons (top). The 
equivalent prompt photon and positron energy spectra in the calorimeter are shown 
in the middle and bottom panels. The low Csl threshold of the trigger discriminators 
is set at about 4.5 MeV. 

Fig. 38. The Netscape experiment control window: different users can remotely mon- 
itor the state of the detector, the status of the data acquisition system, and the cur- 
rent values of the database variables. The collaborators with the password-protected 
access can submit electronic logbook entries, stop or pause an ongoing run or start 
a new run. 

Fig. 39. The radial-azimuthal projection of a /x+ — > e^Vf.v^e^e'' event in the PI- 
BETA single event display (SED, top). The digitized and fitted waveform of the 
incident vr"*" stopping in the active target, followed by the tt"*" IJ-'^t^h decay as seen 
in the SED (bottom). 

Fig. 40. The background-subtracted tt"*" e^Ve energy spectrum in the Csl 

calorimeter. Software gains of all Csl detectors are automatically adjusted by a 
fitting program. The GEANT-simulated detector response is represented by the full 
line histogram. 



62 



Fig. 41. The agreement between the measured /x+ — > e'^i^e'^^ energy spectrum in 
the Csl calorimeter and the Monte Carlo simulation. The data are collected over 10 
runs. The per degree of freedom is 1.3. 

Fig. 42. The calibrated energy spectrum (sum of energies of coincident photon pairs) 
for the final 7r/3 event sample (top). Two-photon energy spectrum produced by Tr~ 
single charge exchange (SCX) on the degrader and active target counters (bottom). 
The rms energy resolution of 4.0 % is consistent with the calorimeter response ob- 
tained for clean 7r/3 events. 

Fig. 43. The timing spread of 220 trigger-defining Csl analog signals at the inputs 
of the UVA 125 discriminator/summer units. The mean time difference between the 
fast proton TDC hit in a Csl module corrected for a time-of-flight delay and the 
TDC value in the reference active degrader counter has rms width of 0.86 ns. 

Fig. 44. Uncorrected (top panel) and linearized (bottom panel) ADC vs TDC scatter 
plots for one representative Csl channel. Displayed data points are accumulated over 
10 consecutive runs. 

Fig. 45. Measured time difference in ns between two coincident photons from 7r° 
produced via SCX in the active target. The top histogram is formed using raw 
TDC values. The bottom histogram demonstrates the effectiveness of the TDC 
linearization. 

Fig. 46. Signal-to-background ratio for the events with one minimum ionizing and 
one neutral particle inside the vr-stop gate. The cuts imposed on the high threshold 
one-arm calorimeter trigger are listed. The top panel corresponds to the vr-stop rate 
of ~ 7.3-10^ 7r+/s while the bottom panel shows the data for the low intensity beam 
of ~ 1.4- lO^TT+Zs. 

Fig. 47. The online time distribution of the tt" events recorded with the 7r/3 trigger 
with respect to the 7r-stop signal. The average value of the two photon TDCs whose 
energies and directions reconstruct to the 7r° is shown. The narrow-peak pattern 
corresponds to the 50.63 MHz radio-frequency cyclotron beam structure. 

Fig. 48. The angular resolution of the PIBETA calorimeter: differences between the 

direction of a minimum ionizing positron calculated using the power algorithm with 
energy depositions in the nearest-neighbor Csl clusters and a reconstructed track 
direction from the MWPC hits. 
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Table 1 

Geometrical specifications of the two cylindrical chambers MWPCi and MWPC2 
used for the charged particle tracking. 



Chamber Parameter 


Inner MWPCi 


Outer MWrC2 


Active length [mm] 


350 


540 


All length [mm] 


580 


730 


Min. diameter [mm] 


90 


206 


Max. diameter [mm] 


152 


255 


Diameter anode [mm] 


120.3 


240.2 


Diameter inner cathode [mm] 


115.2 


235.2 


Diameter outer cathode [mm] 


125.0 


245.2 


Inner cathode-anode gap [mm] 


2.45 


2.5 


Volume 01 chamber [cm ] 


660 


2035 


Total chamber thickness [mg/cm^] 


53.9 


74.8 


iotai chamber thickness [rad. length] 


1.4 • 10-3 


2.0 • 10-3 


Anode wires 






Number of anode wires 


192 


384 


Anode wire spacing [mm] 


1.96 


1.96 


Anode wire tension (min-max) [N] 


0.44-0.48 


0.44-0.51 


Resistance of anode wire [Ohm] 


110 


155 


Resistance between anode-anode [Ohm] 


560k 


560k 


Capacity between anode-ampl. [pFJ 


1000 


1000 


Cathode strips 






Number of inner cathode strips 


2-64 


192 


Angle of slope inner cathode [deg] 


36.960 


44.602 


Number of outer cathode strips 


2-64 


192 


Angle of slope outer cathode [deg] 


-33.650 


-42.339 


Width of cathode strips [mm] 


3.0 


2.4 


Gap between strips [mm] 


0.4 


0.3 


Resistance of strips (6u Al) [Ohm] 


~1.5 


-2.0 
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Table 2 

Optical properties of the pure Csl scintillators used for the PIBETA calorimeter 
(Manufacturers: Bicron Corporation and Harkov Institute for Single Crystals). 



Quantity 


Average Value 


"D ., r 

li(M(M (nice 


Density 


4.53 g/cvo? 


[8] 


Radiation Length 


1.85 cm 


[8] 


Refractive Index at 500 nm 


1.80 


[8] 


Refractive Index at 315 nm 


1.95 


[22] 


Nuclear Interaction Length 


167 g/cTo? 


[8] 


dE/dx\rain 


1.243 MeV/g/cm^ 


[8] 


Photonuclear Absorption (tt^ photons) 


0.675% 


[59] 


Light Attenuation Length 


103 cm 


[58] 


Fast Component Wavelength 


305 nm 


[25,26] 


Fast Component Decay Time 


7 ns 


[25,26] 


Slow Component Wavelength 


450 nm 


[25,26] 


Slow Component Decay Time 


35 ns 


[25,26] 


Lower Wavelength Cutoff 


260 nm 


[22] 


Fast-to- Total Light Output Ratio 


0.76 


[23,31] 


Light Output 


66.2 Photoel./MeV 


[23,31] 


Light Output Nonuniformity 


+0.28 %/cm 


[23,31] 


Light Output Temperature Coefficient 


-1.56 %/°C 


[23] 


Time Rcsohition (wrt PV tag) 


0.68 ns 


[31] 


Radiation Resistance 


See Sec. 18 


[57] 


Stability 


Slightly Hygroscopic 


[25,26] 
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Table 3 

Csl cluster definitions for the PIBETA calorimeter. Sixty overlapping clusters CLOO- 
CL59 each contain nine Csl modules. The first six clusters CL00-CL50 make the 
first Csl supercluster SCOO, and so on. 



C1/CbI 


1 


2 


3 


4 


5 


6 


7 


8 


9 


Cl/Csl 


1 


2 


3 


4 


6 


6 


7 


8 


9 


CLOO 


60 


200 


160 


120 


20 


10 


110 





100 


CLOS 


66 


215 


166 


125 


26 


IB 


115 


5 


105 


CLIO 


70 


120 


170 


180 


130 


30 


20 


154 


80 


CL15 


76 


125 


175 


185 


136 


35 


25 


159 


86 


CL20 


80 


30 


130 


197 


188 


140 


40 


147 


178 


CL25 


85 


35 


135 


192 


183 


145 


45 


142 


173 


CL30 


90 


50 


40 


140 


190 


171 


150 


138 


181 


CP35 


95 


55 


45 


145 


195 


176 


155 


133 


186 


CL40 


100 


110 


10 


60 


150 


161 


210 


121 





CL45 


106 


115 


16 


55 


166 


166 


205 


126 


6 


CL60 





10 


20 


30 


40 


50 


70 


80 


90 


CL55 


5 


15 


26 


35 


46 


55 


76 


85 


96 


CLOl 


61 


201 


161 


121 


21 


11 


111 


1 


101 


CL06 


66 


216 


166 


126 


26 


16 


116 


6 


106 


CLll 


71 


121 


171 


181 


131 


31 


21 


150 


81 


CL16 


76 


126 


176 


186 


136 


36 


26 


155 


86 


CL21 


81 


31 


131 


198 


189 


141 


41 


148 


179 


CL26 


86 


36 


136 


193 


184 


146 


46 


143 


174 


CL31 


91 


51 


41 


141 


191 


172 


151 


139 


182 


CL36 


96 


56 


46 


146 


196 


177 


156 


134 


187 


CL41 


101 


111 


11 


51 


151 


162 


211 


122 


1 


CL46 


106 


116 


16 


56 


156 


167 


206 


127 


6 


CL61 


1 


11 


21 


31 


41 


61 


71 


81 


91 


CL56 


6 


16 


26 


36 


46 


56 


76 


86 


96 


CL02 


62 


202 


162 


122 


22 


12 


112 


2 


102 


CL07 


67 


217 


167 


127 


27 


17 


117 


7 


107 


CL12 


72 


122 


172 


182 


132 


32 


22 


151 


82 


CL17 


77 


127 


177 


187 


137 


37 


27 


166 


87 


CL22 


82 


32 


132 


199 


185 


142 


42 


149 


175 


CL27 


87 


37 


137 


194 


180 


147 


47 


144 


170 


CL32 


92 


52 


42 


142 


192 


173 


152 


135 


183 


CL37 


97 


57 


47 


147 


197 


178 


157 


130 


188 


CL42 


102 


112 


12 


52 


152 


163 


212 


123 


2 


CL47 


107 


117 


17 


57 


157 


168 


207 


128 


7 


CL52 


2 


12 


22 


32 


42 


62 


72 


82 


92 


CL57 


7 


17 


27 


37 


47 


57 


77 


87 


97 


CL03 


63 


203 


163 


123 


23 


13 


113 


3 


103 


CLOS 


68 


218 


168 


128 


28 


18 


118 


8 


108 


CL13 


73 


123 


173 


183 


133 


33 


23 


152 


83 


CL18 


78 


128 


178 


188 


138 


38 


28 


167 


88 


CL23 


83 


33 


133 


195 


186 


143 


43 


145 


176 


CL28 


88 


38 


138 


190 


181 


148 


48 


140 


171 


CL33 


93 


53 


43 


143 


193 


174 


153 


136 


184 


CL38 


98 


58 


48 


148 


198 


179 


158 


131 


189 


CL43 


103 


113 


13 


53 


153 


164 


213 


124 


3 


CL48 


108 


118 


18 


58 


158 


169 


208 


129 


8 


CL53 


3 


13 


23 


33 


43 


63 


73 


83 


93 


CLB8 


8 


18 


28 


38 


48 


58 


78 


88 


98 


CL04 


64 


204 


164 


124 


24 


14 


114 


4 


104 


CL09 


69 


219 


169 


129 


29 


19 


119 


9 


109 


CL14 


74 


124 


174 


184 


134 


34 


24 


153 


84 


CL19 


79 


129 


179 


189 


139 


39 


29 


168 


89 


CL24 


84 


34 


134 


196 


187 


144 


44 


146 


177 


CL29 


89 


39 


139 


191 


182 


149 


49 


141 


172 


CL34 


94 


54 


44 


144 


194 


170 


154 


137 


180 


CL39 


99 


59 


49 


149 


199 


175 


159 


132 


185 


CL44 


104 


114 


14 


54 


154 


KiO 


214 


120 


4 


CL49 


109 


119 


19 


59 


159 


165 


209 


125 


9 


CL54 


4 


14 


24 


34 


44 


54 


74 


84 


94 


CL59 


9 


19 


29 


39 


49 


59 


79 


89 


99 
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Table 4 

Physics triggers built into the PIBETA electronics logic. "PS" labels the prescaled 

trigger events, "D" represents dedicated triggers, "C" stands for the Csl calorimeter, 
"B" is the 7r+ beam. "LT" and "HT" denote a low discriminator threshold {Ejji = 
4.5 MeV) and high discriminator threshold (£^ht = 52.0 MeV), respectively. 



Number/ 


Trigger 


Trigger 


Prescaling 


Typical 


Lal-)ol 




Logic 


Factor at IMUz 


Rate- (llz) 




Two- Arm HT D 


CHgOTrGoB^ 


1.0 


8 




Two- Arm LT PS 


CggOvrGpsoBTT 


512 


6 


3/tL 


One- Arm HT PS 


C^^ovrGpsoB^ 


16 


40 


4/tH 


One- Arm LT PS 


CgOvrGpsoB^ 


32768 


4 


5/Tp 


Prompt HT PS 


Cg o TrSps' 


1024 


2 


6/TL 


Cosmics LT D 


CV o o B^ 


1 


46 




TT-in-Beam HT PS 


timed BCoADoATorf 


1.6 • 10^ 


0.5 


8/TH 


Cosmics LT PS 


CVpsoC^oB^ 


8 


0.5 


9/Te 


e-in-Beam HT PS 


timed BCoADoATorf 


32 


0.2 


10/Tr 


Random Trigger PS 


random 


1 


2 


11/T3S 


Three- Arm LT PS 


CggOvrGpsoBTT 


64 


8 


12/TLH 


Two- Arm LT-HT PS 


C^gLoTrGoB^ 


4 


5 
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Table 5 

Summary of changes in gain factors and energy resolutions of active elements due 
to radiation exposure. 



Detector 


Radiation 


Starting 


Final 


Starting 


Final 




Dose (rads) 


Gain 


Gain 


Res. (%) 


Res. (%) 


BC 


2.0- 


10^ 


1.00 


0.66 


13.1 


13.9 


AD 


1.4- 


10^ 


1.00 


0.75 


7.8 


8.5 


ATo 


5.6 ■ 


10^ 


1.00 


0.76 


7.2 


10.8 


ATi 


6.3 • 


10^ 


1.00 


0.89 


8.4 


8.6 


AT5 


1.5 • 


10^ 


1.00 


0.74 


7.7 


8.9 


PVo 


4.4. 


10^ 


1.00 


0.99 


31.5 


32.1 


PVi 


4.0- 


10^ 


1.00 


0.99 


29.9 


31.1 


PVeto 


4.0 ■ 


10^ 


1.00 


0.95 


26.2 


27.9 


Cslo 


45.9 


1.00 


0.89 


5.1 


5.5 


Csl2 


44.5 


1.00 


0.92 


4.9 


5.2 


Cslii 


116.5 


1.00 


1.02 


6.0 


6.1 


Csli9 


93.0 


1.00 


0.65 


5.0 


5.3 


Csli02 


151.7 


1.00 


0.86 


6.0 


6.5 


Csli65 


125.5 


1.00 


0.74 


5.4 


5.8 


Calo 


118.5 


1.00 


0.83 


5.5 


6.0 



68 




FIGURE 1 



69 




FIGURE 2 



70 




FIGURE 3 



71 




FIGURE 4 



72 



en 


700 


1 1 1 1 1 1 1 1 1 1 1 1 1 


1 1 1 1 1 1 1 1 


c 

D 




p„+=113.4±0.5 MeV/c 




o 
o 


u w w 










• 






500 






CD 
















E 


400 




• 










z 


300 








200 


* 


• 




100 


m 

0' 

^ . • • 




;■. , , 



?12 112.5 113 113.5 114 114.5 115 



7T^ Momentum (MeV/c) 



FIGURE 5 



73 




FIGURE 6 



74 



{ 


• 


» 

10 mm 


®Beam 



FIGURE 7 



75 




FIGURE 8 



76 



400 F 



300 



200 - 



100 








xVndf72'.99 


/ 56 






Constant 


360.4 


± 


4.179 


Mean 


107.1 


± 


0.2743 


Signna 


13.60 


± 


0.1745 







50 100 150 200 250 300 

BC ADC (ch) 



FIGURE 9 



77 



Sideview ^ Frontview 




FIGURE 10 



78 




FIGURE 11 



79 



40 mm 



30 mm 




front view side view 



FIGURE 12 



80 




FIGURE 13 




81 



xVnd034.5 / 207 




To ADC (ch) 

FIGURE 14 



82 




83 




FIGURE 16 



84 



c 



10 - 



O 



0.13% IX 



99.5% 7y' 



< - 



-10 - 



20 



X 



0.37% e' 



□ □ 

□ □ □ □ 



tt""— Stop 
Trigger 



20 



15 



10 



5 10 

rf -ATTDC (ns) 



FIGURE 17 



85 




AzcH (mm) 

FIGURE 18 



86 



75 



50 



25 

1 —Arm Ml e"^ Trocks 





50 100 150 200 250 300 350 

MWPC, (p (deg) 

n ' I ' I ' I ' I ' I ' I ' I : 

75 
50 
25 

1 —Arm Ml e"^ Trocks 

° 50 100 150 200 250 300 350 

MWPC2 (f (deg) 



FIGURE 19 



87 



(i) 




FIGURE 20 



88 



A 



1 —Arm Charged Track 



+ + 



Protons 



V 



1/ 



f 

/ 

>V 







2 



3 



4 5 
PV Energy (MeV) 



FIGURE 21 



89 




FIGURE 22 



90 







<£pv>=99.223±0.002% 



50 100 



l|D 

Entries 

Mean 

RMS 



160d83 
1 15159CH 
180.1 
104.0 



150 200 250 300 350 

(^pv (deg) 



FIGURE 23 



91 




FIGURE 24 



92 



(i) (ii) 






FIGURE 25a 



93 



(vii) 



(viii) 





FIGURE 25b 



94 





FIGURE 26 



95 



btain 



theta 
160 




WO 360 
P™ 



FIGURE 27 



96 



Wavelengthsliifter 



CM CNJ 



20,0 



300,0 



160,0 



2900,0 



3000,0 



steel frame 
steel plate 

corrugated cardboard 
wavelength shifter 
scintillator 




FIGURE 28 



97 



^ 600 F 
o 

< 

o 400 



200 : 













2— Arm HT Trigger 




Runs 21 100^21 109 







'40 -20 20 40 60 80 100 

Linearized CV TDC (ns) 



FIGURE 29 



98 



10000 



7500 



5000 



2500 



Entries 
Mean 



426369 
136.1 



-Arm Cosmic Trigger 
Ep,>10 MeV 




10 



2 



10 



100 



200 



300 400 
Ed (MeV) 



Eci,c2>52 MeV /* \ 



Oc,,c2>150 



Entries 
Mean 



16916 
309.8 





2-Arm CT 
— GEANT 








200 



400 



600 800 
Eci + Ec2 (MeV) 



FIGURE 30 



99 




FIGURE 31 



100 



(U 

> 
o 

00 
Q 

"D 
(D 
N 

E 
o 



1 

0.8 
0.6 

0.4 
0.2 

h- 



































1 




+ 


— in — Degrc 


jder Syster 


n Function 







































































50 



100 



150 



200 250 
Tinne (ns) 



FIGURE 32 



101 




0-' 







Ij 
c. 




CJ- 

'sr 


eoQ 


" O 

Q 
< 


400 


U.. 

to, 


500 




^ICJOO 
"5 

*| 600 

n 

200 



' 1 II I' ■ r- 



^00 



lOOO -2O0 




TOC (n^) 



FIGURE 33 



102 



CO ' ' ' ' ' ' 

c 

i> 10 - * Random Event Trigger 

o *. l„+=800k nVsec 

^ 10^r 



1 n 2 : 



'ti 



10 20 30 40 50 60 

Csl Energy (MeV) 



FIGURE 34 



103 



r timing 



LAM 




Tip Trigger 1 999 with tinree LB 500 



revised 29-Jul-1999 



FIGURE 35 



104 



Tek Sfism sooms/s 



3 Acqs 




A: ISO. On 
@: -56. On 



500mVQ Ch2 SOOmVfi M 25.0ns Chi ^ 
CTT3 1.00 V Ch4 2.00 V 



S40mV 14 Nov 20 
09:30:50 



FIGURE 36 



105 



I 500 

> 
u 

400 

1 300 : 



200 : 

■ + 
100 - 



,1+ 











50 



Entries 


18776 


Mean 


62.23 - 


RMS 


29.84 


Prompt 




Trigger 





ID=Proton 



100 



150 



200 





Entries ' 11746 


_ \ 


Mean 22.57 


RMS 18.05- 


- \ 


Prompt 
Trigger 

ID=Photon 


i + 




i 

* 




+* + 













50 



100 



150 



200 



2000 



1000 



t + 

+ 


Entries ' 40007 
Mean 32.22 
RMS 12.02" 


Prompt 
Trigger 




ID=Positron 


+ 




+ 









50 



100 150 200 

Ec, (MeV) 



FIGURE 37 



106 



File Edit View Go Communicator ' 








j € ^ 3 ^ 

Back Forward Reioad Home 


Search Netscape 


Print Security BMi 


li \^'' Boofcmari<;s ^ Location: ||ittp://H 


.pisi . ch/ 




©1' What's Related 


i; ^ Members ^ Webtvlail ^ Connections 


^ BIzJournal ^ SmartUpdate ^ MMplace | 



IE 



MIDAS ei^erimeiiit "Pibeta" 



Wed Oct 11 20:13:29 2000 



Stop Pause] ODB CHAf| Massages | ELog| Alsiims| Programs | 



Histoiyl Gonfig Help 



Trigger s ettmgs Rates Ratio s Clismb ejs PI D H aiidb o ok WebC am Ac c elerator PAW 



Rtffi #36005 


Rurming 


Alarms: On 


Restart: Yes 


Data dii: /data 


Start: Wed ct 1 1 1 9 :55 :1 7 2000 


Rtinrung time: 0hl8ml2s 


Equ^dnent 


FE Node 


Events 


Event rate[/s] 


Datarate[kB/s] 


Analyzed 


Tngger 


Trigger Frontend@pc812 


36553 


28.3 


53.2 


100.0% 


Scaler 


Trigger Fronteiid@pc812 


108 


0.0 


0.0 


100.0% 


HV 


SC Fronteiid@pc8D9.psi.ch 


19 


0.0 


0.0 


100.0% 


Enviromnent 


3C Fronteiid@pc809.psi.ch 


19 


0.0 


0.0 


100.0% 


Chamber 


Trigger Fronteiid@pc812 


19 


0.0 


0.0 


' 100.0% 


Beamliiie 


3C Fromend@pc809.psi.ch 


2 


0.0 


0.0 


100,0% 


Channel 


Active 


Events 


MByvritten 


GB total 


OrtuiSeOOS.mid 


Yts 


36781 


61 .445 


793.986 


Lazjr Label 


Pcogcess 


File Name 


# Files 


Total 


2000-3S 


100 % 


nm36003.mid 


49 


90.5 K 


Lazjr Destination 


Pcogcess 


File Name 


Speed \khts] 


Total 


psarduve 


100 % 


nm36003.mid 


735.3 


55.0 % 



20:09:0S[Tng£ei: Frontend] [dsc.c:£lS:dsc_spai:cify] Address 21£1, module 2, channel IIS 
belowiange 2170-2179 



Logger [pc2106] 



Anal/ser [po2106] 



La2/_Tape [po2106] 



La2y_FTP [pc2106] 



SC Frontend [pc809.psich] Trigger Frontend [pcSl 2] 



irdittpd [pc2106] 



-3^ -s-ii (Si^ m 



FIGURE 38 



107 



fMercador [^Front [ Side [DSC | 



p ositron 




Tracks 

I3HEB | run3 17 5Q.mid|[NQ: 139Sl| | Run: 3 17"50] | Ser: 13 93 Q|[3^J u I- 2 OOP 1 7:0 7:5 7| 



RBela Single Event Display 



File E^dit C^oTo Jlelp 
fMercadar |'Frant [■'Side I'DSC | 




Channel: \2S3 I 



HHHH | run3 6S70dsc.mi^ | No: 12\ [ Run: 36a70| | Ser: 2 96 | | 2 7- No v- 2 OOP 2 1:3 4:57 | | Trg: OOOOOOOOO 1 OO | 



FIGURE 39 



108 




FIGURE 40 



109 



40000 



^ 35000 



CO 

c 



l_J 



Z 30000 
I 25000 
z 20000 



DATA 
GEANT 




Ee. (MeV) 



FIGURE 41 



110 



in 
c 

>2000 

LJ 



O 



^ 1500 



Entries 


5901 1 


Mean 


127.6 


RMS 


6.401 



1000 - 



500 - 











DATA 
GEANT 



50 




100 



+ V + 

7T T\ e V, 



150 200 
E„o (MeV) 



c 

CD 

^ 6000 - 

O 
i- 

CD 

^4000 - 



2000 



Entries 


198716 


Mean 


128.4 


RMS 


6.786 











sex 7T° AE/E=4.0% 




7i"CH, 



50 



100 



150 200 
E„o (MeV) 



FIGURE 42 



111 



CO 

o 
o 

CD 

CD 
Q 

If) 
O 

CD 
Z5 



30 



20 



10 







Csl Trigger Timing 



0.5 







0.5 




1 1.5 

Csl TDC 



2 2.5 
DEC TDC (ns) 



FIGURE 43 



112 



if) 

^ 10 - 
o 

Q 

— 5 - 
o 

- 



-5 - 

-10- 1 -Arm LT Trigger! Csl No. 33, Uncorrected TDC ^ 
10 20 30 40 50 60 



10 - 

5 - 

- 



-5 



1 - '' "Arm LT Trigger! Csl No. 33, Corrected TDC 

10 20 30 40 50 60 

Csl Energy (MeV) 



FIGURE 44 



113 



ID 





244.5 ' / ' 87 


P1 


45.62 


P2 


1086. 


P3 


-0.2543 


P4 


1.304 




2000 - 

1500- c^=0-84ns 
1000 - 
500 - 




-8 -6 -4 -2 




xVndf' 


199^6 ' / ' 81 


PI 


1.564 


P2 


1899, 


P3 


0.1505E-01 


P4 


0.8414 



2 4 6 8 
TDC 7i - TDC 72 (ns) 



FIGURE 45 



114 



800 



600 



400 - 



Region B 



-e+ 
-C 



>20.0 MeV 



E;>55.6 MeV 



200 




0:,,>4O.O 



xVndf 


l'02!5 ' 7 ' 55 


P1 


607.2 


P2 


-0.5636E-02 


P3 


1.140 


P4 


239.9 


P5 


3.091 


P6 


-0.2432 







10 -7.5 -5 -2.5 



2.5 5 7.5 10 
U (Csl) - t, (Csl) (ns) 



00 - 



50 - 



Region B 
Ee^>20.0 MeV 
E>55.6 MeV 








-10 -7.5 -5 -2.5 



xVndf 


^1.00 ' / '81 


PI 


1 13.2 


P2 


-0.6861E-01 


P3 


0.8837 


P4 


2.028 



2.5 5 7.5 



10 



te. (Csl) - L (Csl) (ns) 



FIGURE 46 



115 



j5 4500 F 900060 

^ - , , Entries 290849 . 

> 4000- ^f=(19-751 ±0.010) ns ^^^^ 

^ RMS 56.28 

O 3500 - 



_g 3000 - 
I 2500 - 
2000 - 
1500 - 
1000 - 
500 - 




-63.10 -23.83 15.75 55.29 94.74 134.13ns 

-83.13 -43.30 ,1 -3.73 35.53 74.93 114.48 



100 -50 50 100 150 

Csl (71+72)72 TDC - DEG TDC (ns) 



FIGURE 47 



116 




2 4 6 8 10 12 



©csi-Qmwpc (deg) 



FIGURE 48 



117 



